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test  data  records  to  generate  the  sixteen  synthetic  maneuvers  used  in  the 
study.  The  synthetic  data  are  corrupted  by  realistic  measurement  noise 
levels  representative  of  the  actual.  The  EBM  method  is  a  two  stage  process 
In  the  first  stage,  the  corrupted  data  are  processed  one  maneuver  at  a  time 
using  an  extended  Kalman-Bucy  filter/ Bryson-Frazier  smoother  to  obtain 
estimated  values  of  the  states  and  of  the  forces  and  moments  acting  on  the 
aircraft.  In  the  second  stage,  subspace  modeling  together  with  a  step-wise 
multiple  linear  regression  technique  is  used  to  process  the  estimated  values 
from  all  sixteen  maneuvers  together  in  parallel  to  identify  a  global  state/ 
control  dependent  model  of  the  three  force  coefficients  and  the  three  momen 
coefficients.  The  identified  global  model  compares  well  with  the  true 
values  of  the  T-2C  wind  tunnel  model.  The  EBM  method  is  demonstrated 
to  accurately  identify  the  high  nonlinearities  of  the  T-2C  wind  tunnel  model. 
The  nrpdiction  aeeuracv  of  the  identified  global  model  is  shown  to  be  excel- 
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Speed  along  aircraft  z-body  axis  with  respect  (m/sec) 

to  a  nonmoving  free  stream  air  mass 

Longitudinal  force  along  x-body  axis  (Newtons) 

Aircraft  c.g.  position  along  the  longitudinal  (m) 

axis  of  earth  fixed  axes  oriented  with  respect 
to  aircraft  body  axes  by  Euler  angles.  Positive 
when  c.g.  is  forward  of  axis  system  origin. 

Distance  between  c.g.  and  sensor  i  measured  (m) 
along  x-body  axis,  i =ar,0,V.  Positive  for  sensor 
forward  of  c.g. 

Distance  between  c.g.  and  longitudinal  accelero-  (m) 
meter  at  i  measured  along  x-body  axis, 
i  =  NOSE.  Positive  for  accelerometer  forward 
of  c.g. 

Distance  between  c. g.  and  lateral  accelerometer  (m) 
at  i  measured  along  x-body  axis,  i  ~  c.g.,  nose, 
tail.  Positive  for  accelerometer  forward  of  c.g. 

Distance  between  c.  g.  and  vertical  accelerometer  (m) 
at  i  measured  along  x-body  axis,  i  -  c.g.  rose 
tail.  Positive  for  accelerometer  forward  of 

c.g. 

Lateral  force  along  y-body  axis  (Newtons) 

Aircraft  c.g.  position  along  the  lateral  axis  of  (m) 
earth  fixed  axes  oriented  with  respect  to  air¬ 
craft  body  axes  oriented  with  respect  to  aircraft 
body  axes  by  Euler  angles.  Positive  when  c.g. 
is  right  of  axis  system  origin. 

Distance  between  c.  g.  and  sensor  i.  Positive  for  (m) 
sensor  right  of  c.  g. 

Distance  between  c.  g.  and  longitudinal  accelerometer  (m) 
at  i  measured  along  y-body  axis,  i  =  NOSE. 

Positive  for  accelerometer  right  of  c.g. 
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Distance  between  c.g.  and  lateral  accelero¬ 
meter  at  i  measured  along  y-body  axis, 
i  =  c.g.,  NOSE,  TAIL.  Positive  for 
accelerometer  right  of  c.g. 

i 

(m) 

yz. 

i 

Distance  between  c.g.  and  vertical  accelero¬ 
meter  at  i  measured  along  y-body  axis, 
i  =  CG,  NOSE,  TAIL,  R.L.  Positive 
for  accelerometer  right  of  c.g. 

(m) 

z 

Vertical  force  along  z-body  axis 

(Newtons) 

1 

Z 

Aircraft  c.g.  position  along  the  vertical  axis 
of  earth  fixed  axes  oriented  with  respect  to 
aircraft  body  axes  by  Euler  angles.  Positive 
when  c.g.  is  below  axis  system  origin. 

(m) 

z. 

l 

Distance  between  c.g.  and  sensor  i  measured 
along  z-body  axis,  i  =  a,  P,  V.  Positive  for 
sensor  below  c.g. 

(m) 

z 

X. 

1 

Distance  between  c.g.  and  longitudinal  accelero¬ 
meter  at  i  measured  along  z-body  axis, 
i  =  NOSE.  Positive  for  accelerometer  below 

c.g. 

(m) 

z 

yi 

Distance  between  c.g.  and  lateral  accelero¬ 
meter  at  i  measured  along  z-body  axis, 
i  =  CG,  NOSE,  TAIL.  Positive  for  accelero¬ 
meter  below  c.g. 

(m) 

z 

z. 

1 

Distance  between  c.g.  and  vertical  accelero¬ 
meter  at  i  measured  along  z-body  axis, 
i  =  c.g.,  NOSE,  TAIL,  R.L.  Positive  for 
accelerometer  below  c.g. 

(m) 
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CL 

Angle  of  attack 

(rad) 

• 

a 

Time  derivative  of  cl 

(rad/ sec) 

Sideslip  angle 

(rad) 

b 

Time  derivative  of  0 

(rad/sec) 

6 

a 

Aileron  deflection  (+  6  -  left  wing  down  roll) 

a 

(rad) 

6 

e 

Elevator  deflection  (+5^—  nose  down  pitch) 

(rad) 

6 

r 

Rudder  deflection  (+  -  nose  left  yaw) 

(rad) 

C'  O 

a  a 

5r’  5r  See  Equations  4. 5. 5  and  4. 5.  6  for 

definition 

9 

Euler  pitch  angle 

(rad) 

P 

Air  density 

(kg/m3) 

<t> 

Euler  bank  angle 

(rad) 

'll 

Euler  yaw  angle 

(rad) 

LIST  OF  SUPERSCRIPTS 

» 

Denotes  intermediate  quantity 

LIST  OF  SUBSCRIPTS 

o 

Denotes  reference  value 

c. 

Center  of  gravity 

L 

Left  wingtip 

m 

Designates  a  measured  quantity  when  used  as  a 
subscript 

R 

Right  wingtip 

s 

Stability  axes 

NOSE 

Nose 

TAIL 


Tail 
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1.  INTRODUCTION 

1.1  BACKGROUND 

The  inability  of  empirical  and  theoretical  techniques  to  accurately  pre¬ 
dict  the  nonlinear  aerodynamics  of  the  high  angle  of  attack /sideslip  domain  for 
aircraft  has  spurred  interest  in  extracting  these  effects  from  actual  flight  data. 
Success  in  this  endeavor  will  allow  predicted  effects  to  be  correlated  with  the 
actual  effects  to  determine  the  scaling,  trending,  reformulation,  etc.  that  must 
be  applied  to  the  predictions  to  match  the  actual.  Performing  this  correlation 
for  several  existing  aircraft  will  build  confidence  in  the  designer's  ability  to 
predict  the  performance  and  stability  and  control  characteristics  of  proposed  air¬ 
craft. 

The  T-2C  is  a  light  jet  trainer  aircraft  which  has  been  instrumented  and 
operated  by  the  Naval  Air  Development  Center  (NADC)  for  the  express  purpose  of 
producing  high  angle  of  attack  flight  data.  The  data  along  with  a  wind  tunnel  model 
of  the  T-2C  has  been  supplied  to  Dynamics  Research  Corporation  (DRC)  for 
processing  with  an  advanced  system  identification  method  called  the  Estimation- 
Before-Modeling  (EBM)  technique.  The  desired  results  of  applying  an  extraction 
method  such  as  the  EBM  technique  to  flight  data  are  (1)  optimal,  smoothed,  time 
history  estimates  of  the  vehicle  state  parameters  and  (2)  state -dependent  models 
for  the  coefficients  of  the  forces  and  moments  acting  on  the  vehicle.  The  bottom 
line  of  such  an  application  is  whether  the  coefficient  models  can  accurately  pre¬ 
dict  vehicle  motion  for  inputs  different  from  those  analyzed  to  develop  the  models. 

The  first  phase  of  the  processing  of  the  T-2C  data  and  wind  tunnel  model 
is  a  simulation  study  to  check  out  the  accuracy  of  the  EBM  technique  using  syn¬ 
thetic  T-2C  data.  Such  synthetic  data  are  generated  by  exciting  the  T-2C  wind 
tunnel  model  with  the  controls  and  initial  states  contained  in  the  actual  T-2C  flight 
data  records.  In  the  second  phase,  the  actual  T-2C  flight  data  are  processed. 

This  volume  details  the  results  of  the  first  phase. 
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1.  2  SCOPE  OF  SIMULATION  STUDY 

This  study  consists  of  the  following  four  parts: 

•  Simulation 

•  Estimation 

•  Modeling 

•  Prediction 

In  the  simulation  part,  sixteen  maneuvers  of  synthetic  T-2C  flight  data 
are  generated  by  exciting  NADC's  T-2C  wind  tunnel  mode?  with  the  control  time 
histories  and  the  initial  states  of  sixteen  actual  T-2C  flight  data  records.  The 
data  are  corrupted  by  the  addition  of  white,  zero -mean  gaussian  noise  having 
realistic  standard  deviations.  In  the  estimation  part  the  corrupted  data  of  the 
sixteen  maneuvers  are  processed  by  a  filter /smoother  algorithm  to  provide  esti¬ 
mates  of  the  states  and  of  the  forces'  and  moments'  time  histories.  In  the  model¬ 
ing  part,  these  estimates  are  processed  by  a  stepwise  multiple  linear  regression 
technique  over  a  network  of  a,  0  subspaces  to  identify  global  state -dependent 
models  for  the  coefficients  of  the  forces  and  moments.  The  identified  models  are 
compared  with  their  counterparts  in  the  wind  tunnel  model.  In  the  prediction  part 
a  new  maneuver  as  well  as  the  original  sixteen  maneuvers  are  used  to  test  the 
accuracy  of  the  identified  model. 

1.3  ESTIMATION-BEFORE -MODELING  (EBM)  APPROACH 

The  Estimation-Before-Modeling  (EBM)  technique  is  a  two-step  approach 
to  system  identification.  The  first  step  is  the  estimation  of  states  and  aerodynam¬ 
ic  forces'  and  moments'  time  histories  using  a  nonlinear  estimator.  In  the  second 
step,  the  state  and  control  input  space  is  divided  into  many  small  subspaces. 
Estimated  states  from  all  maneuvers  that  fall  within  a  chosen  subspace  are  used 
in  modeling  the  force  and  moment  coefficients  over  that  subspace.  For  each  sub¬ 
space  the  state  and  control  dependent  model  is  obtained  using  Stepwise  Multiple 
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Linear  Regression  (SMLR).  The  EBM  approach  provides  the  following  benefits: 

a.  Estimation  in  the  time  domain  provides  better  data  for 
the  modeling  phase. 

b.  No  a  priori  state -dependent  model  need  to  be  selected 
to  perform  state  estimation  and  the  identification  of 
aerodynamic  derivatives. 

c.  Multiple  maneuvers  can  be  processed  in  parallel  to 
identify  aerodynamic  derivatives  with  greater  accuracy. 

d.  The  global  state -dependent  modeling  region  can  be  sub¬ 
divided  into  smaller  regions,  and  the  identification  of 
aerodynamic  derivatives  can  be  performed  within  each 
region  separately. 

1.3.1  Estimation 

The  equations  that  describe  the  motion  of  an  aircraft  can  be  expressed 
as  (i.  e. ,  see  Appendix  A): 

x  =  f  (x)  +  B  (t) 

where  x  represents  the  9-dimensional  state  vector  (u,  v,  w,  p,  q,  r,  9,  0,  z), 

f  (x)  denotes  the  kinematic  terms  and  B  (t)  *  (B ,  (t),  B_(t),  .  .  .  ,  B~  (t),  0,  0,  0) 

with  B^(t)  =  acceleration  in  the  x^- direction  due  to  forces  and  moments  acting  on 

the  aircraft.  The  accelerations  B^^  (t),  i  =  1,  2,  .  .  .  ,  6,  are  linear  functions 

of  the  undimensionalized  aerodynamic  coefficients  C  (t),  C  (t),  C  (t),  C,,(t),  C  (t), 

x  y  z  *  m 

and  C  (t).  These  aerodynamic  coefficients  are  unknown  functions  of  the  states 

n  •  * 

a,  0,  p,  q,  r,  a,  0,  etc.  and  the  controls  6,6,6.  Since  the  models  of 

ear 

C  .....  C  are  unavailable,  the  accelerations  B  (t),  i  -  1,  2,  ....  6,  are 
x  n 

necessarily  unmodeled  accelerations.  The  estimation  problem  of  the  first  step 
is  that  of  estimating  the  state  of  a  nonlinear  dynamical  system  in  the  presence 
of  unmodeled  accelerations. 

One  method  of  attack  on  the  problem  is  the  spline  estimation  method  as 
outlined  in  Ramachandran,  et.  al.  [  l  ]  and  detailed  in  Schneider,  et.  al-  [  2  3- 
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The  approach  taken  herein  is  to  approximate  the  unmodeled  accelerations 
B^t),  i  =  1,  2,  ....  6,  by  Gauss-Markov  processes  and  to  apply  an  extended 
Kalman  filter/ Bryson-Frazier  smoother.  This  method  of  approach  is  described 
in  Appendix  B. 
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1.3.2  MODELING 

The  estimated  values  of  C  ,  C  ,  C  ,  etc.,  along  with  the  states  and 

x  y  z 

control  inputs  from  the  estimation  phase,  are  further  processed  to  remove  the 
thrust  and  c.g.  offset  effects  and  converted  to  suitable  form  for  the  purpose  of 
identification. 


The  steps  involved  in  the  model  identification  are  outlined  below.  A  few 
passes  through  all  these  steps  are  usually  necessary  to  obtain  a  satisfactory  global 
model. 

Construction  of  Subspaces 

The  advantage  of  using  subspaces  for  modeling  is  that  we  can  determine 

the  highly  complex  aircraft  characteristics  accurately,  using  only  simple  models. 

The  first  step  in  the  model  identification  stage  is  the  choice  of  subspaces.  The 

force  and  moment  coefficients  are  generally  functions  of  several  aerodynamic 

variables.  Hence,  it  is  first  necessary  to  determine  which  independent  variables 

will  significantly  affect  the  force  or  moment  coefficient  whose  model  we  want  to 

obtain.  Here  it  is  important  to  bear  in  mind  that  the  objective  is  not  just  to 

arrive  at  a  mathematical  formula  that  will  best  fit  the  data  at  all  points  but  to 

determine  the  aerodynamic  characteristics  of  the  airplane  as  accurately  as 

possible.  The  selection  of  significant  independent  variables  and  the  nature  of  their 

influence  on  the  dependent  variable  will  determine  the  choice  of  subregions. 

For  example,  from  previous  experience,  available  data  on  the  airplane  and 

engineering  judgement  it  is  known  that  for  the  given  airplane  Cm  is  influenced  very 

strongly  by  a  and  6  ,  less  strongly  by  0,  by  q  in  a  linear  fashion  and  very  weakly 

by  lateral  controls.  In  this  case  the  subspaces  will  be  bound  by  z,  5  and  f},  with 

e 

small  widths  for  z  and  6g  and  larger  widths  for  0.  Having  settled  upon  the  importart 
variables  and  the  choice  of  variables  to  define  the  boundaries  for  modeling  subspaces, 
the  boundaries  for  these  regions  can  be  chosen  by  the  nature  of  the  phenomenon 
being  studied  and  known  aircraft  characteristics.  For  instance,  if  we  are  analyzing 
the  aircraft  behavior  in  the  post  stall  regime,  we  may  like  to  use  very  narrow 
regions  for  u  right  after  stall  and  larger  regions  farther  away  from  stall  angle. 
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Also  it  is  very  important  that  the  width  of  the  subregion  in  each  independent  variable 
should  be  much  greater  than  the  standard  deviations  of  the  estimation  error  and 
measurement  noise  for  that  independent  variable. 

Having  roughly  determined  the  subspaces,  the  next  step  is  to  judge  the 
adequacy  of  the  data  and  the  data  density  distribution  in  each  region.  This  is 
done  by  simply  determining  the  total  number  of  data  points  that  fall  within  each 
subspace  using  the  data  from  all  maneuvers.  Only  the  maneuvers  that  give 
consistent  and  high  quality_results  in  the  estimation  stage  are  used  in  all  stages  of 
identification.  A  comparison  of  the  total  number  of  points  in  each  subspace  and 
the  number  of  independent  variables  to  be  used  in  the  model  will  give  us  an  idea 
as  to  the  sufficiency  or  scarcity  of  data  for  obtaining  a  model  in  each  region.  This 
will  also  help  us  to  ascertain  whether  the  subspaces  and  the  type  of  model  will 
have  to  be  redefined.  A  high  density  may  make  it  possible  to  use  still  narrower 
regions  to  obtain  more  accurate  models  with  lower  order  models.  Further 
refinements,  where  necessary,  could  be  made  based  upon  the  quality  of  the 
identified  model  and  the  statistics  of  the  independent  variables  in  each  sub¬ 
space. 

Selection  of  Independent  Variables 

The  estimated  state  and  control  variable  time  histories  obtained  by  the 
stochastic  estimation  method  are  used  to  compute  the  independent  variables 
for  the  regression.  It  is  possible  to  choose  any  desired  number  of  independent 
variables  as  candidates  for  the  model.  Each  independent  variable  is  formed  by 
combining  one  or  more  of  the  following  terms  or  higher  powers  of  them: 


(a  -  aQ)  ,  <0  -  0Q)  ,  (se  -  5e  )  ,  (6a  -  5a  )  ,  (ir  -  5r  )  ,  , 


rb  o_c  £b_ 

2 V  '  2V  '  2 V  * 


Engineering  judgment,  known  aerodynamic  properties  of  the  aircraft 
and  the  nature  of  the  problem  under  investigation  can  be  used  to  eliminate  super- 
flous  terms. 


llT-fi-J 
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Model  Determination 

To  model  a  given  coefficient  its  value  along  with  the  appropriate  states 
and  control  inputs  are  chosen  from  all  maneuvers  if  the  corresponding  a, 

/3  and  control  surface  deflections  fall  within  the  chosen  subspace.  By  using 
information  from  all  maneuvers  it  can  be  ensured  that  maximum  information  is 
available  for  modeling. 

The  SMLR  technique  is  employed  in  each  subspace  for  several  values  of 
F-critical  ratio  and  a  model  is  determined  for  each  F-critical  ratio.  The 
adequacy  of  a  model  is  determined  by:  a)  the  mean  and  standard  deviation  of  the 
dependent  and  indoDendent  variables  in  the  subspace,  b)  the  noise  in  the  measure¬ 
ments  ,  c)  the  statistics  on  the  estimates  of  the  states  and  forces  from  the 
stochastic  filter,  d)  the  F-ratio  and  standard  deviation  of  each  term  of  the 
determined  model,  e)  the  Residual  Sum  of  Squares  (RSS)  and  multiple  correlation 
coefficient  R  ,  and  f)'  engineering  judgment. 

The  best  model  for  a  given  region  is  obtained  by  applying  the  preceding 
considerations  and  analyzing  the  models  obtained  at  each  step  of  the  regression 
for  each  F-critical  value. 


The  aerodynamic  parameters  determined  from  all  subspaces  are  assembled 
to  provide  the  global  model. 
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SECTION  2  j 

|  T-2C  WIND  TUNNEL  MODEL  -AND  SYNTHETIC  T-2C  D-AT-A 

:  For  the  purpose  of  checking  out  the  two  step  approach,  synthetic  data  are 

|  generated  from  the  wind  tunnel  model.  The  wind  tunnel  model  was  supplied  by  j 

j  NADC  [7].  Details  of  the  wind  tunnel  model  and  the  computation  of  the  total 

I  aerodynamic  force  and  moment  coefficients  along  the  aircraft  body  axes  are  j 

!  i 

|  given  in  Appendices  A  and  D.  The  normal,  axial  and  side  forces  are  given  by 


X  = 

PV2SC  +  T  COS  ; 

2  x  -e 

(2.1) 

Y  = 

T pv2scy 

(2.2) 

Z  = 

(2.3) 

The  angle  l ^  between  the  X  force  body  axis  and  the  thrust  axis  was  supplied 
as  zero.  The  roll,  pitch  and  yaw  moments  are  given  by: 


_  P'/2Sb 


(2.4) 


M 


N 


p  V2Sc 

-  -T^cm  ♦iis-T»8y-il(ZtTSmy 

2  n 


+  T  cos  £  l 
e  z 


(2.5) 


(2.6) 


L 2=1. 
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In  order  to  help  make  the  simulation  be  more  like  the  real  flight  data,  ; 

the  mass,  control  inputs,  trim  velocity  and  trim  altitude  are  taken  to  be  the 

same  as  those  used  in  the  actual  flight  test  of  the  airplane.  The  geometry  of 

the  aircraft,  mass  and  inertia  characteristics  were  furnished  by  N.ADC  and  are  , 

given  in  .Appendix  A.  .All  the  aircraft  constants  are  the  same  for  all  maneuvers 

except  the  mass  which  is  given  for  each  individual  maneuver.  No  thrust  history  ! 

was  available  nor  the  value  of  the  distance  between  the  aerodynamic  center  and  | 

the  center  of  gravity.  In  view  of  possible  discrepancies  between  the  wind  tunnel  ; 

model  and  the  actual  aircraft  it  was  necessary  to  calculate  the  trim  angle  of 

attack  for  the  simulation.  For  each  maneuver  the  trim  angle  of  attack  a  .  , 

&  trim 

trim  thrust  T  and  the  distance  between  the  center  of  gravity  and  aerody¬ 

namic  center  2^  were  calculated  to  satisfy  the  trim  equations: 


u  =  v  =  $=  ,S  =  p  =  q  =  r  =  p  =  q  =  r  =  i:,5  =  *  =0  (2.7) 

3l  IT 


Vim  ‘  2 trim 


(2.8) 


Assumption  of  perfect  trim  and  lateral  symmetry  of  the  given  model 
were  used  to  write  the  above  equations. 

The  trim  angle-of-attack  to  balance  the  forces  in  the  vertical  direction 
can  be  calculated  by  solving 


'h  =  qu  -  pv  +  g  cos  6  cos  <+  +  —  =  0 


(2.9) 


and  equation  (2.8)  in  an  iterative  fashion. 
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Having  computed 

compute  T.  . 

trim 


a,  .  .  the  equation  for 

trim 


u 


given  below,  is  solved  to 


u  =  rv  -  qw  -  g sin Q  +  X/m.  =  0  (2.10) 

With  a..  and  T.  .  ,  f  is  computed  by  solving  for 

trim  trim  x  c  J  6 

k  =  T^[M  +  pr(I  -  I  )  +  (r2  -  p2)I  ]  =  0  (  2.11) 

1  Z  X  aZ 

y 


An  inspection  of  the  flight  data  revealed  that  for  several  maneuvers  the 
aileron  and  rudder  deflections  were  not  zero  but  of  small  magnitude  at  trim. 
These  values  were  treated  as  biases  and  subtracted  from  the  flight  time 
histories  to  give  control  inputs  for  the  simulation. 

As  pointed  out  above  the  simulation  checkout  was  to  resemble  as  much  as 
possible  the  system  identification  of  an  actual  aircraft.  For  this  reason  the 
three  controls  £>a>  £>g  and  together  with  the  initial  conditions  of  the  states 
a,  j9,  V,  p,  q,  r,  6,  <t>  and  z  were  copied  from  actual  flight  data  records  and 
used  in  the  generation  of  the  synthetic  data.  Sixteen  maneuvers  of  available  T-2C 
actual  flight  test  data  were  used  for  this  purpose.  These  are  summarized  in 
Table  2-1.  The  first  letter  S  of  SF1M1  denotes  synthetic,  the  F  denotes  flight 
number  and  the  M  denotes  the  maneuver  number  of  that  flight.  .  The  designation 
SFIM1  denotes  that  the  controls  and  the  initial  conditions  used  to  generate  this 
synthetic  maneuver  were  copied  from  maneuver  1  of  flight  1  of  the  actual  flight 
test  data  provided  by  NADC.  The  type  of  controls  is  described  in  the  fourth 
column  for  each  maneuver.  Plots  of  the  controls  for  all  maneuvers  are  given 
in  Appendix  D. 
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Table  2-1  MANEUVERS  USED  IN  THE  STUDY 
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The  number  after  SF  indicates  NADC  flight  number  and  the 
number  after  M  indicates  NADC  maneuver  number. 
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With  the  given  aerodynamic  model  and  the  initial  conditions  calculated  by 
the  above  mentioned  procedure,  the  equations  of  motion  were  integrated  while 
exciting  the  airplane  with  control  inputs  from  flight  test  data  for  the  chosen 
maneuver.  When  this  was  done,  it  was  found  that  a  few  maneuvers  diverged 
after  a  few  seconds,  and  the  responses  for  most  maneuvers  were  totally  different 
from  actual  flight  data.  The  wind  tunnel  model  exhibited  frequencies  and 
amplitudes  that  were  very  much  higher  than  flight  data.  The  plots  within  Figure 
2.  1  give  a  comparison  between  the  response  of  the  wind  tunnel  model  and  the 
aircraft  for  maneuver  15.  It  may  be  noticed  that  the  aircraft  motion  did  not  go 
beyond  an  angle  of  attack  of  2  5  degrees,  exhibited  small  oscillations  about  the 
maximum  and  returned  to  trim  value  when  the  elevator  was  returned  to  trim 
setting,  whereas  the  angle  of  attack  for  the  wind  tunnel  model  first  reached  a 
peak  of  22  degrees  and  continued  to  grow  with  the  amplitude  reaching  nearly  45 
degrees.  Even  after  the  return  of  elevator  angle  the  growth  continued  with  the 
angle  of  attack  going  up  to  60  degrees.  An  inspection  of  the  lateral  states  reveals 
much  larger  magnitudes  and  a  much  greater  degree  of  instability  for  the  wind 
tunnel  model. 

Subsequently  by  a  trial  and  error  approach  the  aerodynamic  tables  were 
modified  slightly  to  provide  more  reasonable  responses.  The  modifications 
consisted  of  very  small  adjustments  in  the  tables  in  the  stall  and  the  post  stall 
regions  and  fixing  the  value  of  Cn  at  -0.06.  The  constant  value  of  -.06  for 
Cn  was  found  to  give  stability  to^the  simulated  data  and  to  provide  data  which 
agreed  better  with  the  actual  flight  data.  For  example  note  the  improvement  in 
the  response  of  the  modified  wind  tunnel  model  of  maneuver  15  as  shown  in 
Figure  2.2.  For  the  modified  case  the  state  a  and  $  return  to 

trim  values  when  the  control  inputs  return  to  trim  values.  The  rates  p,  q  and 
r  become  small  at  trim  control  for  the  modified  case  whereas  they  continue  to 
grow  with  large  amplitudes  for  the  original  model.  The  maximum  values  of 
a,  j3,  p,  q,  r,  9  and  4>  are  40°,  20°,  120  deg/sec.  50  deg/sec,  30  deg/sec,  80° 
and  700°,  respectively  whereas  for  the  original  model  they  are  60°,  50°, 

320  deg/sec,  80  deg/sec,  100  deg/sec,  160°  and  1300°,  respectively.  The 
modified  T-2C  aerodynamic  model  was  used  as  the  basis  for  the  simulation  study 
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Response  of  Modified  Model  for  Maneuver  15 
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and  is  referred  to  as  the  T-2C  wind  tunnel  model  throughout  the  report. 

The  nonlinearity  exhibited  in  the  original  is  similar  to  the  nonlinearities 

contained  in  the  other  dynamic  derivatives  (e.g.  °C,  ,  C  ,  and  C  ).  Plots  of 

*-  **  n 

p  r  r 

and  C  ^  are  given  in  Figures  4.11  and  4. 12.  A  plot  of  C  is  given  in 
p  r  r 

Figure  4. 18.  Note  that  the  modified  wind  tunnel  value  of  is  plotted  as  the 
constant  -.06  in  Figure  4.17.  Consequently,  fixing  the  valui  of  C  was  no  loss  ! 

^  I 

to  checking  out  the  EBM  system  identification  technique  since  there  was  a  gain 
in  stability  and  since  the  other  dynamic  derivatives  provided  similar  nonlinearities J 
for  testing  the  goodness  of  the  EBM  technique.  | 


I 
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3.  FILTER/SMOOTHER  ESTIMATION  RESULTS 

For  each  of  the  sixteen  simulation  maneuvers  discussed  in  Section  2 
the  time  history  of  the  states  were  fed  into  the  measurement  equations  of 
Appendix  A  to  produce  noise-free  measurements.  These  measurements  were 
corrupted  with  white,  zero-mean  gaussian  noise  having  the  standard  deviations 
as  supplied  by  NADC.  These  standard  deviations,  which  are  shown 
in  Table  3-1  in  the  column  under  the  heading  "Measurement  Noise  Level",  were 
derived  from  analysis  of  actual  T-2C  flight  data.  The  corrupted  measurements 
were  processed  using  the  filter /smoother  algorithm  described  in  Appendix  B. 

Table  3-1  also  contains  a  column  for  the  RMS  values  of  the  estimation 
errors.  These  RMS  values  are  obtained  by  taking  the  difference  between  the 
noise -free  measurements  and  the  measurements  obtained  through  substituting 
the  estimated  states  into  the  measurement  equations.  The  RMS  values  in  Table 
3-1  are  for  Maneuver  4;  the  results  for  the  other  maneuvers  are  given  in  Tables 
E-l  through  E-3  of  Appendix  E. 

The  RMS  values  of  the  estimated  errors  of  the  accelerations  (t), 
i  =  1,  2,  ....  6,  are  presented  in  Table  3-2  for  Maneuver  4  and  in  the  Tables 
E-4  and  E-5  of  Appendix  E  for  the  other  maneuvers.  Time  histories  of  true 
versus  estimated  states  and  acceleration  are  plotted  in  Figure  3-1  through  3-7. 
The  true  values  are  plotted  with  a  solid  line  while  the  estimates  are  plotted  with 
a  dashed  line  containing  one  cross  symbol  per  second. 


Estimated 

Quantity 


Units 


RMS  Value  of 
Estimation  Error 
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Using  Table  3-1  we  make  the  following  observations  of  the  estimation 
errors  of  Maneuver  4.  The  filter /smoother  algorithm  has  reduced  the  error 
level  from  .  24  to  .  12  for  the  acceleration  measurement  A  X  NOSE.  The  errors 
on  the  acceleration  measurements  in  the  z -direction  (i.e.,  AZCG,  ....  AZL) 
have  been  reduced  by  a  factor  of  3  or  4.  The  errors  on  the  acceleration  measure 
ment  in  the  y-direction  have  been  reduced  by  a  factor  of  3.  The  errors  on  9 
and  0  have  been  reduced  by  a  factor  of  3.  The  errors  on  the  roll,  pitch  and 
yaw  rates  have  been  reduced  by  less  than  a  factor  of  2.  The  errors  on  air  speed 
are  less  by  a  factor  of  7.  The  errors  on  a  and  0  have  been  cut  down  by  a 
factor  of  5  or  6.  These  observations  hold  for  all  maneuvers  with  the  exception 
of  Maneuvers  1  and  7.  The  RMS  value  of  the  estimation  error  on  the  roll  rate  is 
slightly  higher  than  the  noise  level  for  these  two  maneuvers  in  Tables  E-l  and  E-2 

Making  use  of  Table  3-2  and  Figures  3-5  through  3-7,  we  observe  that 
the  errors  in  the  estimates  of  EL  (t),  i  =  1,  2,  .  .  .  ,  6,  are  about  one-thirtieth 
of  the  peak  values  of  EL  (t).  Moreover,  the  estimate  of  EL  (t)  followed  accu¬ 
rately  the  time  history  of  EL  (t). 


o 

o 


T I  ME ( SEC  ) 

Figure  3-1.  Comparison  Between  the  True  and  Estimated  Values  of 
a  and  {3  of  Maneuver  4 


3-5 


TRUE  VALUE 


I 


* 


o 

o 


Figure  3-3.  Comparison  Between  the  True  and  Estimated  Values  of 
Q  and  R  of  Maneuver  4 
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ESTIMATED  VALUE 


5. OS  82.53 


90.02  97. SO  104.98 
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112.47  119.95 


TRUE  VALUE 


:  .05  82. S3  90.02  97.50  104.98  112.47  119.95 

T I  ME ( SEC  ) 


Figure  3-4.  Comparison  Between  the  True  and  Estimated  Values  of  9  and  0 
of  Maneuver  4 
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SECTION  4  j 

I 

MODEL  IDENTIFICATION  AND  COMPARISON  OF  RESULTS  ! 


Model  comparisons  of  force  and  moment  coefficients  can  be  done  in  either 
j  equation  form  or  in  the  look  up  tables  form.  Both  these  methods  are  equivalent 
j  and  equally  good.  To  illustrate  both  the  approaches,  the  tabular  form  was 
j  chosen  for  longitudinal  models  and  the  equation  form  for  lateral  models. 

i  4.1  CONSTRUCTION  OF  SUBSPACES  FOR  LONGITUDINAL  IDENTIFICATION 

The  method  and  factors  governing  the  choice  of  subspaces  were  discussed 
|  in  Section  1.3.  The  stall  angle  of  attack  for  the  aircraft  was  around  14°.  Since 


Cx>  Cz  and  C^  are  strongly  dependent  upon  a,  |j3|and  5^  the  subspaces  were 
defined  by  the  following  subdivisions: 


5-  -25°.  -15°.  -10°.  0°.  10°  and  15° 


8  :  0°,  2°,  5°.  10°.  15°,  20°  and  25° 


a  :  -8°,  -4°.  0°.  4°.  8°.  12°.  14°.  15°.  16°.  18°.  20°.  25°.  30°,  35°. 
40°  and  45° 


I  The  symmetry  of  longitudinal  coefficients  with  respect  to  sideslip  angle  was 

I  taken  advantage  of  in  using  the  same  model  for  both  positive  and  negative  sideslip 
j  angles.  For  convenience  and  ease  of  comparison  against  the  wind  tunnel  model, 

|  the  subdivision  values  of  a,  [j3|  and  Se  were  chosen  to  be  the  points  where  the 
!  values  of  CT  .  C„  and  C  were  specified  for  the  wind  tunnel  model.  For  real 
|  flight  data  a  few  trials  may  be  necessary  before  settling  upon  the  best  choice 
;  based  on  the  identified  models  in  these  trials.  Also  the  assumption  of  symmetry 


of  C  ,  C  and  C  with  respect  to  3  is  only  a  convenience  to  reduce  the  number 
x  z  m 

of  subspaces  and  increase  the  data  for  modeling  in  each  subspace.  When  asymmetr; 

is  suspected  the  procedure  can  be  used  in  exactly  the  same  way  by  modeling  in  the 
positive  and  negative  sideslip  angle  regions  separately. 


I 
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Using  all  the  maneuvers,  data  points  for  all  subspaees  were  calculated. 

Table  4-la  shows  the  data  distribution  in  the  region  contained  by  ol  =  -4°,  ?/5°, 

!fl  |  =  0°,  30°  and  6  =  -15°,  -10°.  Such  a  table  immediately  shows  us  how 

e 

the  data  points  are  distributed  and  in  what  regions  additional  data  must  be 

generated  for  complete  modeling.  Similar  tables  were  generated  for  all  other 

subdivisions  of  6  and  the  results  have  been  summarized  in  Table  4- lb.  Table 
e 

4-lb  lists  modeling  regions  containing  subspaces  with  at  least  15  points.  The 
number  of  subspaces  for  each  region  is  also  given.  Each  row  of  the  table  shows 
the  maximum  and  minimum  on  6g,  |  £  j  and  a  within  which  all  the  points  were 
located  and  the  resultant  number  of  subspaces.  For  example,  the  seventh  row 
of  the  table  shows  that  for  -15°  <  6^  s  -10°  and  0°  £  |/3  ]  <2°  the  data  points 
were  located  for  angles  of  attack  between  0  and  30°.  As  we  can  see  from 
Table  4-la  this  resulted  in  ten  subspaces  since  we  have  ten  subdivisios  of 
angle  of  attack  for  ol  between  0°  and  30°.  It  was  noticed  that  there  was  only  a 
scattering  of  points  outside  the  regions  given  in  Table  4-lb.  In  general,  it  was 
found  that  towards  the  edge  of  the  regions  the  data  points  were  low  and  for  5° 
s  1/3  I  s  10°  the  data  were  concentrated  below  \j3  |  =  8°. 

From  the  Table  4-lb  we  see  that  the  subspaces  containing  enough  data 
for  modeling  may  not  form  a  covering  of  the  entire  space  and  there  may  be 
holes  (e.g.  third,  fourth  and  fifth  rows  of  the  table).  We  also  see  that  the 
controls  used  were  inadequate  to  generate  data  in  several  regions,  especially 
at  high  j9  and  also  for  the  combination  of  low  ol  and  high  6  .  We  also  do  not 
have  enough  data  points  for  angle-of-attack  beyond  35  degrees. 


cm 
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Table  4-lb  DLT-4  DISTRIBUTION  FOR  LONGITUDES'LL 

MODELING 


'e  . 
nun 

(deg) 

s 

'e 

max 

(deg) 

3 

min 

(deg) 

0 

max 

(deg) 

a 

min 

(deg) 

a 

max 

(deg) 

No.  of 
Subspaces' 

-10 

0 

0 

2 

-4 

25 

10 

2 

5 

0 

25 

9 

5 

10 

4 

14 

3 

5 

10 

16 

18 

1 

5 

10 

20 

25 

1 

10 

15 

8 

12 

1 

-15 

-10 

0 

2 

0 

30 

10  > 

2 

5 

8 

30 

8  i 

5 

10 

14 

30 

6  ! 
i 

-17 

-13 

10 

15 

20 

25 

i 

i 

-25 

-15 

0 

2 

8 

35 

9 

2 

5 

12 

35 

8 

5 

10 

12 

14 

1 

5 

10 

16 

30 

4 

-25 

-30 

0 

2 

14 

35 

7 

2 

5 

16 

35 

5 

5 

10 

18 

35 

4 
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!  From  Table  4-lb  we  see  that  we  obtain  a  total  of  88  subspaces  for 

i 

1  modeling.  Instead  of  frying  to  match  the  entire  data  with  one  equation 

containing  a  very  large  number  of  terms«we  have  divided  the  data  into  88  groups 
such  that  we  can  fit  a  simple  curve  for  each  region  and  combine  these  88  simple 


'  local  models  to  determine  the  global  model. 


As  described  in  Section  1.3  ,  the  SMLR  method  chooses  the  best  model 
to  fit  the  data  based  on  the  supplied  value  of  F  -  Critical  Ratio.  The  computer  ! 
program  of  Reference  [4]  was  suitably  modified  for  aircraft  parameter  identifi-  1 

cation.  Since  the  longitudinal  coefficients  depend  mainly  on  a  ,  8,  q,  6  and 

•  .  e.  ! 

a.  as  a  first  try  it  was  decided  to  use  linear  models  in  0£,  p,  0  ,  q  and  a  in  all  j 

e 

regions  and  try  more  complex  models  if  necessary.  I 


j  In  each  subspace  the  reference  point  is  fixed  at  the  subspace  corner 

corresponding  to  a  .  ,  1 3  !  .  and  5  and  q  =  a  =0.  The  coefficients  C  . 

1  r  0  min  ^  min  e  x 

max 

i  C  and  C  are  then  expressed  as  a  first  order  Taylor  series  expansion  about  the 
;  z  m 

!  reference  point.  The  equation  for  C  can  be  expressed  as 
!  x 

'  C  =  C  +  c  (a  -  a  )  +  C  (  1/3  !  -  \8  *  )  +  C  (5  -  5  ) 

i  x  Xo  xa  0  x8  0  X5  e  eo 

1  e 


C  +  C 

x  2\  x: 

q  a 


(4.2.1) 


!  where  a  ,  3  !  and  c  are  values  of  a  ,  \8  !  and  6  at  the  reference  point  and  , 

o  o  e  e  , 

|o 

1  C  is  the  value  of  C  at  the  reference  point.  Similar  equations  can  be  written  , 
x  x  1 

o  ; 

for  C  and  C  . 

i  z  m  i 

:  | 

1  Identification  was  performed  in  all  subspaces  and  the  models  were 

I  analyzed.  Generally  except  for  C^,  the  linear  equations  were  adequate  to 


follow  the  variation  of  force  and  moment  coefficients  with  high  multiple 
correlation  coefficients.  An  inspection  of  the  individual  terms  showed  that 
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for  most  regions  the  standard  deviations  of  the  parameters  of  the  model 
were  low.  However  in  the  regions  near  the  boundaries  where  there  were 
less  points  the  uncertainties  in  the  estimates  of  the  parameters  were  somewhat 
high  and  in  some  cases  the  multiple  correlation  coefficient  was  also  low.  It  was  J 
noticed  that  the  identified  model  contained  terms  of  Cm^  or  Cm&  ,  but  their 
values  were  nearly  the  same.  When  both  Cm^  and  Cm£.  were  forced  into  the 
model  (corresponds  to  F-critical  values  being  zero  )  the  estimated  value  of  Cm 


and  C  were  far  different  from  their  values  obtained  for  other  F-ratios  but 
m  . 
a 


their  sum  was  equal  to  the  value  of  Cm  +  C 

mq  m 


a 


This  indicated  that  C„  and 

mq 


Cm  ,  cannot  be  obtained  independently  and  only  their  sum  can  be  estimated. 


a 

The  same  was  observed  for  CZq  and  Cz^.  also.  Comparison  of  individual  terms 
of  the  models  for  all  regions  showed  rapid  changes  in  CZ(^  and  Cm^  after  stall. 


For  C  and  C  the  nonlinear  term  q  a  was  included  and  a  was  dropped.  It 
z  m 

was  also  found  that  because  of  larger  noise  on  unrealistic  values  were 

obtained  for  C  .  Consequently  only  a,  8  and  6  were  used  for  modeling  C  . 

x  e  x 

q 

Identification  was  repeated  in  all  regions.  Models  were  obtained  with 

high  multiple  correlation  coefficients  and  low  residuals,  except  for  C  .  The 

2  x 
value  of  R  was  low  for  C^  even  though  the  RSS  was  within  the  noise  level  of 

stochastic  filter  estimates.  The  multiple  correlation  coefficients  for  C z  were 

above  0.  9  for  most  cases  and  between  0.  96  and  0.  99  for  C  .  Further  addition 

m 

of  higher  orders  terms  in  a  and  8  and  cross  product  of  and  a  made  no 
further  improvement  indicating  the  adequacy  of  the  model.  In  a  few  cases  when 
the  number  of  data  points  was  large,  the  multiple  correlation  coefficient  was 
increased  by  reducing  the  6^  width  and  splitting  the  individual  subspace  into  two 
supspaces.  The  correlation  coefficients  for  the  Cx  model  were  in  the 
interval  of  0.7  -  0.8  or  sometimes  even  lower.  It  was  found  that  the 


standard  deviation  of  the  error  on  the  estimate  of  Cx  from  the  estimation  proeess  j 
was  nearly  of  the  same  order  as  the  standard  deviation  of  Cx  itself  in  those 
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regions.  (This  is  due  to  a  higher  noise-to-signal  ratio  in  the  measurement 

of  C^).  Using  higher  order  terms  in  the  models  did  not  change  the  model  at 

all  indicating  that  the  low  correlation  was  more  due  to  noise  than  the  structure 

of  the  model  itself.  As  we  will  see  later  in  Section  4.3,  the  excellent  agreement 

between  the  estimates  and  wind  tunnel  values  for  all  these  coefficients  in  almost 

all  the  regions  also  pointed  to  this.  To  test  the  sensitivity  of  the  identified  model 

to  spurious  terms,  the  lateral  states,  p,  r,  5  and  5  were  included  as  addition- 

a  r 

al  candidates.  Identification  was  again  performed  for  C  ,  C z  and  Cm  in  many 

regions.  It  was  found  that  the  regression  algorithm  successfully  rejected  the 

extraneous  terms  in  almost  all  the  cases.  However,  when  the  multiple  correlation 

coefficient  for  Cx  was  very  low,  one  or  more  of  the  lateral  states  were  also 

included  in  the  C  model. 

x 

Table  4-2  gives  an  example  of  the  coefficients  of  Cm 

model  for  6^  between  -15°  and  -10°  and  jjS'between  5°  and  10°.  As  pointed 

out  earlier  the  data  points  fell  within  angles  of  attack  of  14  and  30  degrees. 

Hence  no  models  exist  for  a  greater  than  30°  or  a  less  than  14°.  We  see  that  noi 

only  the  values  of  the  coefficient  Cm  but  also  its  derivatives  agree  with  the 

wind  tunnel  model.  We  have  essentially  chosen  the  regions  to  be  small 

enough  such  that  the  nonlinear  behavior  of  the  force  and  moment  coefficients 

can  be  adequately  and  accurately  represented  by  simpler  linear  models.  As 

such.the  derivatives  C  ,  C  ,  C  etc.  represent  the  best  average  values 

m  m.  m„ 

°  6e  5 

of  these  derivatives  over  the  given  subspace.  Hence  there  may  be  slight 

discontinuities  in  the  values  of  C  ,  C  and  C  on  the  subspace  boundaries 

x  c  m 

(each  subspace  has  six  boundary  surfaces).  Similar  tables  were  obtained  for 
all  coefficients  in  all  regions  and  these  tables  were  used  to  build  a  global 
model, 

4.3  LONGITUDINAL  GLOBAL  MODEL 

Having  constructed  models  for  individual  regions  our  next  task  was  to 
assemble  these  to  form  a  global  model  which  can  be  used  for  comparison 
against  the  wind  tunnel  model  and  to  predict  the  response  of  the  estimated  model. 
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Table  4-2  Example  of  Estimated  Coefficients  for  C 

Model  Equation 


(— !  5  <  6  ^  _ 10 


5°  ^  |/3[s  io°) 


a 

(deg) 


14- 15 

15- 16 

16- 18 
18-20 
20-25 
25-30 


c 

m 

o 

C 

ma 

C 

m6e 

Wind 

Wind 

Wind 

‘  ! 

Estimated 

Tunnel 

Estimated 

Tunnel 

Estimated 

Tunnel 

-0.0397 

-0.0311 

-1.0000 

-1.0439 

-1.7660 

-1.710 

a 

i 

-0.0463 

-0.0451 

-1.3369 

-1.6492 

-1. 5775 

-1.680 

1 

i 

-0.0749 

-0.0710 

-1.5224 

-1.6519 

-1.6204 

-1.640 

-0. 1246 

-0.1210 

-1. 531 

-1.3464 

-1.4352 

-1.460 

1 

. 

i 

-0.1460 

-0.1690 

-1.758 

-1.5864 

-0.9715 

-1.047 

l 

-0.2829 

-0.2872 

-0.7236 

-0. 5999 

-0.7444 

-0.782 

1 

• 

\ 

C 

m/3 

C 

mq 

C 

m 

qa 

a 

(deg) 

Wind 

Wind 

Wind 

Estimated 

Tunnel 

Estimated 

Tunnel 

Estimated 

Tunnel 

14-15 

-0.0880 

-0.1287 

-14.98 

-13.20 

315. 15 

431.22 

15-16 

-0.1658 

-0.1552 

-8.59 

-7.70 

624. 57 

676.26 

16-18 

-0.1566 

"0. 1956 

1.82 

3.20 

22.92 

132.31 

18-20 

-0.0502 

-0. 1433 

3.60 

2.41 

-1.26 

-15.09 

20-25 

-0.0320 

-0.U456 

1.66 

2.3 

-47.39 

-42. 184 

25-30 

0.0 

-0.0089 

-2.92 

-1.835 

-27.85 

-70.88 

& 
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It  may  be  recalled  that  in  general  even  though  data  points  lay  between  0°  and  35° 

of  angle  of  attack,  the  concentration  of  these  points  depended  on  angle  of  sideslip 

and  elevator  deflection.  For  instance  for  6  between  -15°  and  -10°  and  B 

e 

between  5U  and  10u  we  could  obtain  models  ?nly,  a  between  14  and  35  degrees. 
Table  4-lb  gave  the  regions  for  which  data  points  were  available  for  model¬ 
ing. 

Models  from  more  than  one  adjacent  subspace  can  be  used  to  fill  the  hole 
in  the  table.  Hence,  engineering  judgment  and  the  quality  and  behavior  of  the 
model  in  each  subspace  was  taken  into  consideration  in  projecting  the  values  of 
force  and  moment  coefficients  into  regions  where  no  models  were  available  for 
lack  of  sufficient  data.  With  this  approach  C  ,  C z  and  Cm  values  were  extrapola¬ 
ted  at  various  points  to  complete  the  tables  for  0°  s  a5-  35°,  0°  =»  1/3 1  s  10°  and 

-25°  £  6  £  0°.  It  should  however  be  remembered  that  caution  must  be  exercised 
e 

in  using  the  extrapolated  values.  The  reason  is  that  such  an  extrapolation  may 
not  be  accurate  if  the  variation  of  the  force  or  moment  coefficient  is  highly 
nonlinear  and  the  models  vary  drastically  from  one  subspace  to  the  next. 
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4.4  DISCUSSION  OF  LONGITUDINAL  MODELING  RESULTS 

The  longitudinal  modeling  results  are  discussed  in  this  section.  Identified  j 

values  of  C  ,  C z  and  are  compared  against  wind  tunnel  values  in  Figure  4.1, 

4.2  and  4.3  respectively.  The  values  of  C  ,  C  and  C  are  plotted  as  functions  of 

x  z  m 

angle  of  attack  for  elevator  settings  of  0°,  -10°,  -15°  and  -25°  at  0  =  0°,  2°,  5° 

and  10°.  The  dynamic  derivatives  C  +  C  and  C  +  C  are  compared 

z  z  •  m  m  * 

q  a  q  a 

against  wind  tunnel  values  in  Figures  4.4  and  4.5  respectively.  Excellent  results 
were  obtained  by  the  identification  technique.  As  previously  mentioned.  Appendix 
D  presents  a  brief  description  of  the  wind  tunnel  model. 

Figure  4. 1  shows  a  comparison  of  estimated  Cx  for  various  values  of  /J 

and  6  for  the  identified  model  and  the  wind  tunnel  model.  At  0  =  0°  the  C 
e  x 

increases  as  the  angle  of  attack  increases  reaching  a  maximum  at  about  14-15 
degrees.  Then  there  is  a  very  sharp  decline  in  the  value  due  to  stall.  This  is 
followed  by  a  gradual  increase  or  decrease  depending  on  the  elevator  angle  and 
angle  of  attack.  We  see  that  the  identification  technique  has  been  successful 
in  matching  these  variations  extremely  well.  We  see  similar  trends  for  C x  at 
other  j3  also.  We  also  observe  very  close  agreement  between  the  estimates  and 
wind  tunnel  values  for  j3  =  2°  and  0  =  5°  for  all  a  and  6^.  For  0  =  10°  also  we 
get  excellent  fit  on  Cx  except  at  ot  =  15°.  This  was  due  to  lack  of  data  points 
around  0  =  10°,  a  =  15°  and  extrapolation  from  neighboring  subspaces. 

i 

| 

Figure  4.2  gives  a  comparison  of  C  for  the  identified  model  and  the  i 

z 

wind  tunnel  model  for  various  a,  0  and  &  .  in  general,  with  an  increase  in 

e  I 

angle-of-attack  the  normal  force  coefficient  increases  (negatively)  in  a  near  ' 

o 

linear  fashion  up  to  a  =  12  .  It  continues  to  increase(negatively)  up  to  stall  where 
it  drops  and  continues  to  decrease  up  to  18  or  20  degrees.  Then  it  starts 
increasing  (negatively)  again  though  with  a  much  lower  slope  than  in  the  ; 

pre-stall  region.  The  actual  values  and  slopes  vary  depending  on  0  and 
We  find  excellent  agreement  between  the  identified  model  and  the  wind  tunnel 
model  for  all  a,  0  and  &  . 
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The  Ditching  moment  coefficient  C  provides  a  much  stronger  test  for  I 

m  1 

the  identification  method  in  view  of  its  greater  sensitivity  to  o  .  Figure 

4.3  compares  C  for  identified  and  wind  tunnel  models  at  various  a,  3  and  i  . 

m  e 

An  inspection  of  these  figures  show  that  not  only  but  also  the  control  effectiver 

ness  C  and  the  pitch  moment  slope  C  depend  on  a,  3  and  5  . 

m*  c  m.  r  e  : 

-  a 

e 


In  general  Cm  varies  linearly  up  to  an  angle  of  attack  of  12  degrees  after 
which  it  decreases  more  rapidly  up  to  25°  degrees.  Beyond  2  5°  there  is  very  j 

j 

little  change  in  Cm  with  change  in  cl.  The  pitch  moment  slope  is  nearly  zero 
between  25  and  30  degrees  for  several  cases.  The  control  effectiveness  shows 
a  great  degree  of  nonlinearity.  At  low  angles  of  attack,  the  control  effective¬ 
ness  becomes  more  negative.  The  control  effectiveness  between  -10°  and  -15°  ; 

is  only  a  fraction  of  that  between  0°  and  -10°.  Beyond  -15°  there  is  no  control  ! 
effectiveness  at  all.  With  increase  in  angle  of  attack  the  control  effectiveness 

between  6  =0°  and  6  =10°  and  also  between  -  =  -10°  and  5  =  -15°  decreases.  I 

e  e  e  e 

However,  for  elevator  angle  below  -15°,  the  control  effectiveness  is  negligible  ! 

at  low  angle  of  attack.  The  actual  values  vary  depending  on  a,  3  and  5  .  ; 

e 


The  dynamic  derivatives  C  +  C  and  C  +  C  did  not  vary 

z  z .  mm-  J 

q  a  q  a 

significantly  with  3  and  for  the  identified  models.  Also  the  given  wind 

tunnel  models  assumed  no  variation  of  these  derivatives  with  respect  to  3  and 

5  .  Hence  for  each  a  the  average  values  of  C  +  C  and  C  +  C  were 
e  z  z-  m  m* 

q  a  q  a 

computed  using  the  estimated  values  from  all  subspaces.  Figure  4.4  presents 

a  comparison  of  estimated  C  +  C ^  .  for  the  identified  and  wind  tunnel  models. 

q  a 

The  value  of  C  +■  C  stays  nearly  constant  at  -6.0  up  to  angle  of  attack  of  14°. 

zq  za 

At  stall  the  derivative  jumps  to  a  positive  value  of  1. 55  between  a  =  14  and 

cl  =  16°.  Beyond  a  =  16°  it  gradually  decreases  to  -2.  Even  though  influence 

of  C  on  C  was  minimal,  the  identified  model  followed  the  wind  tunnel  values 
z  z 

q 

quite  well  except  at  a  =  14°. 


4-19 


«cc 


DYNAMICS  RESEARCH  CORPORATION 


SYSTEMS  DIVISION 


4.5  PRELIMINARIES  FOR  LATERAL  MODELING 


In  Section  1.3  the  considerations  for  choosing  subregions  for  modeling 

were  discussed.  For  the  lateral  case  the  coefficients  are  influenced  by  both 

aileron  and  rudder  deflections.  It  may  be  advantageous  to  divide  the  deflection 

of  these  surfaces  into  several  subdivisions.  But  this  can  lead  to  an  extremely 

large  number  of  subspaces  for  modeling  ,  since  the  total  number  of  subspaces 

will  be  given  by  N  x  Nfl  x  N.  x  N,  where  N  ,  N0,  N.  and  N.  are  the 
apoo  a  pc  6 

a  r  a  r 

numbers  of  a,  ft,  6  and  6  divisions  respectively.  To  take  advantage  of  the 

3.  r 

antisymmetric  nature  of  lateral  forces  and  increase  the  data  density  in  subspaces, 
divisions  of  control  inputs  were  done  indirectly  by  defining  a  new  set  of  control 
variables.  Hence  it  was  necessary  to  use  only  i  and  f$  for  defining  regions  for 
modeling.  To  simplify  the  modeling  procedure  it  was  also  decided  to  assume  no 
interference  between  the  effects  of  aileron,  rudder  and  elevator  and  that  the 
effects  of  aileron  and  rudder  can  be  directly  added. 

For  a  given  angle  of  attack,  because  of  the  symmetric  configuration  of  the 
aircraft,  we  can  assume  the  lateral  force  coefficients  to  be  antisymmetric  with 
respect  to  the  sideslip  angle  when  the  control  deflections  are  zero.  Stated  form¬ 
ally. 


C  (  -/3)  -  -  C  (fi ) 

y  y 


(4.5.1) 


CJL  (-P)  =  -  Ct  «*> 


(4.5.2) 


C  (  -0)  =  -  C  03) 
a  n 


(4.5.3) 


The  static  force  coefficient  at  any  a,  0,  6  and  6  can  be  viewed  to  be  gener- 
ated  by  first  placing  the  aircraft  in  the  corresponding  angles  of  attack  and  side¬ 
slip  with  respect  to  the  wind  and  then  deflecting  the  controls  to  their  respective 
values;  that  is. 


C  =  C  (a,  0)  +  A  C  (6  )  +  A  C  (6  ) 
y  y  y  a  y  r 


(4.5.4) 


Similar  equations  can  be  written  for  C .  and  C  . 

*  n 
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At  a  given  angle  of  attack,  if  we  use  rudder  deflection  alone  the  following 
four  configurations  are  possible: 

(i)  0  positive  and  6^  positive 

(ii)  0  positive  and  6  ^  negative 

(iii)  0  negative  and  5^  positive 

(iv)  0  negative  and  6^  negative 

Now  let  us  consider  the  changes  produced  in  the  lateral  forces  using  the  values 

at  0  ~  0°,  6  =  0°  and  6  =  0s  as  the  reference  values.  Because  of  the  anti- 

a  r 

symmetric  nature  of  lateral  forces,  the  changes  in  the  force  coefficient  due  to 
configuration  (i)  will  be  equal  in  magnitude  and  opposite  in  sign  to  the  changes 
produced  by  (iv).  Similarly  the  change  due  to  configuration  (ii)  will  be  equal  in 
magnitude  and  opposite  in  sign  to  that  of  (iii).  We  notice  that  the  magnitude  of 
contribution  is  the  same  if  the  sign  of  the  product  0  6^  is  the  same  irrespective 
of  the  individual  signs  of  0  and  6  .  This  is  true  for  the  effect  of  aileron  also. 


With  the  above  considerations  in  mind  let  us  define  new  control  variables 


6=6  ,  0  6  2  0 
r  r  r 


5  =  0  .  0  8  *  0 

r  r 


(4.5.5) 


5=0,  06  <  0 

r  r 


6"  =  6  ,  0  6  <  0 
r  r  r 


Similarly, 
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=  6a  * 

=  0  , 

=  0  , 

>  5a<0 

6a* 

*Sa<0 

(4.5.6) 


We  see  that  6^  corresponds  to  positive  when  0  is  positive  and  to  negative 
6^  when  0  is  negative  and  that  one  configuration  is  the  mirror  image  of  the  othei] 
This  is  true  for  the  configurations  corresponding  to  6^,  6*  and  also 

Using  Equations  (4.5. 1)  through  (4.  5.  6)  we  can  write  Cy  as 

Cy  *  Cy  (a,  J^|)  +  A  Cy  (6+)  +  A  Cy  (6;>  +  A  Cy  (6+) 

+  A  Cy  (§a)  +  A  Cy  (Dynamic)  ,  0  z  0 
when  A  Cy  (Dynamic)  is  the  contribution  of  dynamic  derivatives. 

Let  us  assume  that  6  ,6  ,8  ,6  remain  the  same  and  0  just  re- 

verses  the  sign.  Recalling  the  definitions  of  6*  ,  8^,  6*  and  6^  and  the  anti¬ 
symmetric  nature  of  the  lateral  forces  we  can  write 


(4.5.7) 


Or, 


C  =-C  (a,  ]j3 1 )  -  A  C  (6+)  -  A  C  (6~)  -AC  (6+) 
y  y  yr  yr  ya 

-  A  Cy  (6  & )  +  A  Cy  (dynamic) 


-C  =C  (a,  I/JD  +  a  C  (6+)  +  AC  (O+iC  (8+) 
y  y  yr  yr  y  a 

+  A  C  (6  )  -  A  C  (dynamic),  0  <  0 
y  a  y 


(4.5.8) 


(4.5.9) 
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Similar  equations  can  be  written  for  C ,  and  C  •  From  the  above 

*  n 

equations  we  see  that  whenever  0  is  negative  if  we  reverse  the  signs  of  C  ,  C 

y  * 

and  and  the  corresponding  contributions  from  dynamic  derivatives  we  can  map 
these  in  the  positive  0  domain  and  increase  the  data  density  for  modeling. 

It  is  important  to  note  that  the  assumption  of  antisymmetry  was  made  only 

to  facilitate  modeling  by  increasing  the  data  density  in  a  subspace.  Such  an 

assumption  can  be  discarded  when  necessary  and  modeling  can  be  done  separately 

in  the  positive  and  negative  0  domains.  In  such  a  case  there  is  no  need  to  define 
.  + 

6  ,  etc. 

r 

4.  6  IDENTIFICATION  OF  LATERAL  MODELS 

As  noted  in  the  preceding  section,  the  subregions  for  lateral  modeling 
were  created  by  using  a  and  0.  Since  no  control  variable  was  used  directly 
to  define  the  subspaces,  higher  data  density  was  obtained  than  for  the 
longitudinal  case.  This  made  it  possible  to  use  narrower  widths  in  u.  The  sub¬ 
spaces  were  created  by  partitioning  a  at  0°,  4°,  6°,  8°,  10°,  12°,  13°,  14°,  15°, 

16°,  17°,  18°,  19°,  20°,  22.5°,  25°,  27.5°,  30°  and  35°.  Since  few  points 
were  found  when  the  magnitude  of  sideslip  angle  exceeded  ten  degrees  the  divi¬ 
sions  for  |0 1  were  at  0*,  2®,  5°  and  10®.  There  were  also  data  points  in  the 
regions  given  by:  (a)  -4°  s;  u<  0°  and  0°  s  |0|  <2°;  (b)  8°  s  as  12°  and  10°  s: 

|0|  s  15°  and  (c)  20®  <;  a  s  25°  and  10°  s  |0|  £  15®.  As  mentioned  in  Section  4.1 
the  data  points  were  fewer  at  higher  angles  of  attack  and  also  most  of  the  points 
in  the  regions  given  by  5®  s  |0|  s  10®  fell  below  |0  |  =8®.  There  were  not 
enough  data  points  between  a  =  0®  and  a  =  6®  for  |0|  >  5®.  The  total  number  of 
subregions  used  to  obtain  lateral  model  were  55.  Instead  of  trying  to  find  one  very 
complex  model  for  the  lateral  forces,  we  have  reduced  it  to  finding  55  simpler 
models  in  various  regions  and  compiling  them  to  give  us  a  global  model. 

4  -  +  -  * 

As  a  first  try  linear  models  in  a,  0,  6  ,  6  ,  6  ,  5  ,  p,  r  and  0  were  co  i- 

3  3  r  r 

sidered.  From  Eqs.  (4.5.  7)  and  (4.5.  9)  we  see  that  the  data  points  in  the  negative  0 
region  can  be  considered  along  with  the  data  points  in  the  positive  0  region  if  we 
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reverse  the  sign  of  and  the  dynamic  terms.  Then  the  regression  equation 
can  be  written  as 


C=C  +  C  ,0*0 

y  ystat  ydyn 


(4.6.1) 


-  C  =  C  -  C  ,  (3  <  0 

y  ystat  ydyn 

where  C  and  C  are  given  by 


(4.6.2) 


C  =C+C  (a  -  a  )  +  C  (|0|  -  []3  \h-  c  ^  6+ 

\r  Tr  1  1  1  1  tt  4-  n 


ystat  yo  ya 


° 


°  y6  +  a 


+C  6+C.6  +C  6 

V  3  .  V  r  yt'  r 

a  r  r 

C  =  (C  p  +  C  r  -  C  i) 
v  v  v  ys  2V 


ydyn 


(4.6.3) 


(4.6.4) 


The  reference  values  aQ  &  I0Q|  are  the  values  of  u  and  |i3|  at  the  lowest 
comer  of  the  modeling  subregion.  The  regression  equations  for  and 
are  similar.  The  SMLR  method  chooses  the  best  values  of  the  coefficients  in 
the  regression  equations  to  match  C  ,  etc.  in  each  region.  There  may  be  small 
discontinuities  in  etc.  at  the  region  boundaries  since  continuity  is  not  impose* 
at  the  boundary.  However  the  differences  will  be  negligible  if  small  region  widths 
are  used,  the  noise  on  the  stochastic  estimates  of  the  forces  is  small  and  the 
model  matches  the  data  very  closely. 

Identification  was  performed  to  obtain  models  in  all  regions.  The  proce¬ 
dure  was  the  same  as  that  used  for  the  longitudinal  model  determination.  The 
models  seemed  adequate  and  matched  the  forces  with  very  high  correlations  for 
Cjt  and  Cn,  but  not  for  C^.  It  was  also  found  that  the  dynamic  derivatives  changed 
rapidly  after  stall.  Thus  the  nonlinear  terms  pa  and  ra  were  chosen  to  be 
additional  states  for  modeling. 
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The  new  expression  for  Equation  (4.  6.4)  then  becomes 


C  -  (C  P  +  C  r  +  C  p(a-a0)+<^rjr(a-<t0)  ~Cy.fi) ^  (4. 6.  5) 

ydyn  yp  ^r  ^pa  p 


Equations  for  C,  and  C  are  similar.  Using  the  new  model  structure 
dyn  ndyn 

models  were  once  again  obtained  fcrC  ,  and  C  in  all  regions.  Models  were 
obtained  with  high  R  ,  low  RSS  and  low  standard  deviations  for  parameter  values. 
However  in  the  outer  regions  where  there  were  less  points  the  uncertainties  in 
the  estimates  of  the  parameters  were  somewhat  high  and  in  some  eases  the  multi¬ 
ple  correlation  coefficient  was  also  low. 

For  C  (even  though  the  RSS  was  within  the  error  level  of  the  stochastic 

y  2  2 

estimate  of  C  the  R  value  was  lower.  In  general,  for  C,  the  R  values 

y  2  _  2 

were  close  to  one.  For  C  the  R  values  varied  from  0.  96  to  0.  99.  The  R 

n 

values  on  C  however  varied  from  0.  7  to  0.9.  As  was  in  the  case  of 

y 

C  the  reason  once  again  was  that  the  error  on  C  estimates  was  comparable  to 
x  y 

the  variation  of  in  the  region.  (This  is  due  to  higher  noise-to-signal  ratio  in 
the  measurement  of  a  .)  Using  higher  order  terms  did  not  change  the  model  at 

y 

all  indicating  that  the  low  correlation  was  more  due  to  noise  than  the  structure  of 
the  model  itself.  Comparisons  against  the  wind  tunnel  value  verified 
this.  As  in  the  case  of  longitudinal  models, extraneous  terms  such  as  q  and 
were  correctly  discarded  when  introduced  as  additional  variables  for  lateral 
modeling. 
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4.  7  GLOBAL  MODEL  AND  DISCUSSION  OF  LATERAL  RESULTS 

A  short  discussion  of  the  lateral  modeling  results  is  given  below.  In  view 
of  the  assumed  antisymmetry  of  lateral  forces,  the  results  are  presented  for  posi¬ 
tive  or  absolute  values  of  sideslip  angle  only. 

Identified  values  of  C,C,C  ,C  ,C  ,  C  ,  C,  ,  C  ,  C  , 

y  1  *6  lb~  p  r  n  n3 

a  a  r 

C  ,  ,  C  ,  C  and  C  are  plotted  as  functions  of  i  for  3  =  0°,  2°,  5°  and 
n.  +  n  -  n  n 

o  6  p  r 

r  r  r 

10°  and  compared  against  wind  tunnel  values  in  Figures  4-6  through  4-18.  Identi¬ 
fied  values  of  C  ,C  ,,  C  ,  C  ,  ,  C  ,  C  ,  C .  .  ,  C .  ,  C.  ,  C  „ 

V  y<-  y- 


r  6 


C  ,  C  ,  ,  C  ,  C  and  C  are  compared  against  wind  tunnel  values  for 
i  ns+  n.-  n  n  r  6 

n  6  6a  ra 

a  a 

various  a  and  j3  in  Tables  F-2  through  P'-16.  Very  accurate  lateral  models 
were  obtained  using  the  two-step  approach. 


Figure  4. 6  shows  a  comparison  between  the  estimated  and  wind  tunnel 
values  of  C^.  The  values  at  /3  *  10°  were  computed  using  the  model  for  the 
region  5°  s  |0|  £  10°.  With  increasing  0  we  see  a  higher  degree  of  nonlinearity 
in  the  stall  region.  The  identified  model  follows  the  wind  tunnel  values  very  well 
throughout.  In  this  and  all  the  following  curves  for  lateral  coefficients,  due  to 
lack  of  data  for  modeling,  the  values  at  a  =  0°  and  a  =  4°  were  extrapolated 
from  neighboring  regions  when  5°  s  j/3  (  s  10°.  It  was  found  that  for  the 

model  the  most  influential  derivative  was  C  whose  values  are  presented  and 

y|3 

compared  against  the  wind  tunnel  values  in  Table  F-  2.  The  general  trend  is  that 

C  becomes  less  negative  as  angle  of  attack  increases,  reaches  a  maximum 
yjS 

value  around  stall  angle  and  then  becomes  more  negative.  The  value  of  C 


depends  on  both  a  and  0 .  Good  estimates  were  obtained  for  C 

It  was  found  that  the  rudder  deflection  had  a  weak  effect  on  C  .  In  several 

y 

regions  the  contribution  from  6^  was  not  strong  enough  for  lack  of  excitation  of 
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Figure  4.  6  Comparison  of  the  Identified  Model  with  Wind  Tunnel  Values  for  C 
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c+  and  5  or  the  effect  was  not  stronger  than  the  noise  level  on  C  .  The 
r  r  y 

C  ,  and  C  derivatives  and  their  wind  tunnel  values  are  given  in  Tables  F-3 

y‘+  V 

r  "r 

and  F-4.  In  all  tables  a  value  of  zero  is  given  for  a  parameter  when  its  contribu¬ 
tion  was  considered  insignificant  by  the  SMLR  algorithm.  The  derivatives 
C  ,  .  C  and  C  had  much  weaker  influence  on  C  and  values  for  these 

y6  y5-  yr  y 

a  '"a 

could  be  extracted  only  in  certain  regions.  These  are  given  in  Tables  F-5,  F-6 
and  F-7. 


In  Figure  4.  7  we  find  that  the  estimated  values  of  C ^  follow  the  wind 
tunnel  values  very  well.  In  general  with  the  increase  in  angle  of  attack 
increases  in  magnitude  reaching  a  maximum  in  the  stall  region.  With  further 
increase  in  angle  of  attack  it  decreases  in  magnitude.  The  shape  of  the  C 
curve  varies  with  £.  Figure  4.  8  gives  a  comparison  between  estimates  of  C 

0 

and  wind  tunnel  values.  The  parameter  C,  varies  in  a  very  complex  fashion  with 

0 

z  and  #3 . 


Decreasing  aileron  effectiveness  in  roll  is  noticed  with  increasing  angle 
of  attack  in  Figures  4.  9  and  4. 10.  Since  the  values  are  constant  in  a  modeling 
region  they  have  been  plotted  at  the  mid  point  of  the  regions.  Much  of  the  loss 
in  control  effectiveness  occurs  in  the  stall  region.  At  higher  angles  of  attack 
the  control  is  still  effective  though  the  values  are  very  small  compared  to  those 
at  low  angles  of  attack.  Excellent  estimates  have  been  obtained  for  both 
C  and  C.  .  It  was  also  found  that  C.  and  C.  were  not  equal  at 

^  "  g"*"  6~- 

'a  va  a  a 

higher  values  of  a  and  j3.  From  Tables  F-8  and  F-9  we  see  that  compared 

-.he  aileron,  the  rudder  has  an  almost  insignificant  influence  on  C^.  In  some 
-  “.-.o'-.s  the  effect  of  rudder  was  negligible  compared  to  other  terms  in  the  model 
•  's'lmates  were  obtained  for  C.  ,  and  C.  .  In  general  the  estimates  of 

~r  ~ 

5  0 

r  r 

■»r; vatives  were  quite  satisfactory  considering  their  weak  effect. 
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The  roll  damping  derivative  stays  constant  up  to  an  angle  of  attack 

P 

of  ten  degrees  and  starts  decreasing  after  that  (Fig.  4. 11).  Immediately  after 
stall  it  jumps  to  an  unstable  value  of  0.  3  and  stays  positive  up  to  an  a  of  20°. 
Beyond  that  the  roll  damping  is  very  small.  The  identified  model  follows  this 
highly  nonlinear  variation  extremely  well.  As  the  value  of  C  did  not  vary 

Xt 

p 

significantly  with  0,  the  plot  gives  an  average  value. 

The  variation  of  C ^  is  also  very  nonlinear  with  a  with  most  of  the 

x* 

changes  occurring  between  angles  of  attack  of  14°  and  22.5°  (Fig.  4. 12).  Once 

again  we  have  an  excellent  estimate  of  the  coefficient.  The  plot  gives  an  average 

value  of  C  as  there  was  no  significant  change  in  C  with  0.  Even  though  0 
*P  r 

was  used  as  an  additional  variable,  the  SMLR  correctly  excluded  the  0  term 

from  the  identified  models.  The  wind  tunnel  value  for  C.  was  zero.  For  the 

le 

given  model  the  identification  technique  was  successful  in  separating  the  parame¬ 
ters  C  and  C  .  The  estimated  values  of  C  ,  C  and  C  are  given 

Xjq  Hi,  Z 

r  p  a  pa  ra 

in  Tables  F-10,  F-ll  and  F-12. 


Figure  4. 13  compares  estimated  and  wind  tunnel  values  of  Cn«  For 
0  =  2°  and  0  =  5°,  initially  decreases  in  magnitude  with  a,  increases 
near  stall  angle,  decreases  rapidly  after  stall  becoming  nearly  zero  at  a  =  20° 
and  then  remains  very  small  in  magnitude.  For  0  =  10°,  decreases  in 
magnitude  with  w,  falls  sharply  in  the  stall  region  and  reaches  a  negative  value 
at  a  =  18°.  With  further  increase  in  a,  the  value  of  Cn  increases  slightly. 
Very  good  agreement  is  noticed  between  the  estimates  and  the  wind  tunnel  values. 
From  Figure  4 . 14  we  see  that  the  modeling  technique  obtained  good  estimates  of 
C  despite  its  complex  variation.  The  parameter  C  is  positive  at  low  angles 
demines  in  value  rapidly  through  stall  and  becomes  extremely  small  or  negative 
in  the  post  stall  region. 
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Figures  4. 15  anH  4. 16  show  rudder  effectiveness  as  a  function  of  angle 
of  attack  at  various  sideslip  angles.  In  general  there  is  a  decrease  in  rudder 
effectiveness  with  increasing  angle  of  attack.  Much  of  the  loss  in  control  effec¬ 
tiveness  occurs  in  the  stall  regions.  Very  good  estimates  have  been  obtained 
for  both  C  and  C  .  The  influence  of  aileron  was  very  small  compared 

V  V 

r  r 

to  that  of  rudder.  Because  of  this,  only  in  a  few  regions  the  parameters  Cn  + 

6^ 

and  C  were  significant  enough  to  be  estimated  (Tables  F-13  and  F-14). 


From  Figures  4. 17  and  4. 18  we  see  that  the  estimated  values  of  the 
dynamic  derivatives  C  and  C  are  very  close  to  the  wind  tunnel  values. 

As  there  was  little  variation  of  C  and  C  with  /3  the  plots  show  average 

np  nr 

values  of  these  two  parameters.  The  parameter  C  decreases  with  increasing 

r 

angle  of  attack,  experiences  sudden  loss  after  stall  and  becomes  positive.  With 

further  increase  in  a ,  the  parameter  becomes  negative  again  though  the  value 

is  much  smaller  than  at  low  angles  of  attack.  Just  as  in  the  case  of  C.  the 

P 

identification  technique  was  successful  in  separating  C  and  C  #  and  correctlj 

n  Ho 

r  p 

determined  the  parameter  C  .  to  be  insignificant.  The  wind  tunnel  value  for 

0 

C  was  given  as  zero.  Values  of  C  and  C  are  given  in  Tables  F-15  and 
n/j  na  nra 

F-16.  The  parameter  C  was  estimated  to  be  zero  or  negligible  in  all 

nP(r 

regions.  Models  for  the  regions  not  covered  by  Figures  4.  6  through  4. 18  and 
Tables  F-2  through  F-16  are  given  in  Table  F-17. 
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5.  PREDICTION  OF  RESPONSES  FOR  PROCESSED  MANEUVERS 

Having  obtained  the  global  model,  our  next  objective  is  to  assess  how 

well  the  identified  model  predicts  the  maneuvers  used  in  modeling.  The  6 

DOF  equations  of  motion  (Appendix  A)  are  integrated  using  the  control  inputs 

from  each  maneuver.  The  longitudinal  coefficients  C  ,  C  and  C  are 

x  z  m 

calculated  by  interpolation  from  Figures  4. 1  through  4.3  and  adding  the  contri¬ 
bution  of  damping  derivatives  given  in  Figures  4.4  and  4.5.  Equations  (4.5.5), 
(4.5.6),  (4.6.1),  (4.6.2),  (4.  6.  3)  and  (4.  6. 5)  are  used  to  calculate  C  ,  C. 

y  1 

and  C^.  The  parameters  for  these  equations  are  obtained  from  Figures  4. 6 
through  4. 18  and  Tables  F-2  through  F-17.  All  other  parameters  (C  etc. )  were 
negligible.  All  the  equations  of  motion  are  integrated  simultaneously  to  i 

calculate  the  resnonse.  When  a  trajectory  goes  into  a  region  where  the  model 
is  unknown,  estimated  values  ot  states  are  substituted  until  the  trajectory 
returns  to  the  region  where  the  models  are  available. 

Excellent  agreement  between  true  response  and  the  response  predicted 
by  the  identified  model  is  obtained  for  most  of  the  maneuvers.  Figure  5. 1 
shows  the  fit  obtained  for  Maneuver  4.  Comparison  of  responses  for  all  other 
maneuvers  are  given  in  Appendix  G. 

A  comparison  of  predicted  responses  against  true  responses  for  Maneuver;; 
1  through  3,  5  through  10,  first  half  of  Maneuver  11  and  Maneuvers  12  through 
16  are  given  in  the  plots  of  Appendix  G.  The  few  discontinuities  in  the 
states  in  some  of  the  figures  are  caused  by  the  switch  to  the  estimated  values  of 
the  states  in  regions  where  no  models  are  available. 

The  time  histories  of  forces  for  the  predicted  response  also  followed 
those  of  true  responses  very  well.  Figure  5.  2  shows  a  representative  plot.  The 
plot  is  from  Maneuver  4, 

For  several  maneuvers  the  trajectory  went  into  regions  where  no  models 
are  available.  Extrapolating  the  identified  model  into  these  regions  led  to  a 
response  that  is  entirely  different  from  the  true  response.  Hence  estimated 
values  of  states  arc  used  in  these  regions. 
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DYNAMICS  RESEARCH  CORPORATION 

SYSTEMS  DIVISION 

Because  of  the  unstable  nature  of  the  model  in  the  post  stall  region  the  re¬ 
sponses  are  very  sensitive  to  small  deviations  in  the  states.  An  example  of  this  sensi¬ 
tivity  is  shown  in  Appendix  G,  Fig.  G.  1  for  the  second  half  of  Maneuver  11. 

Figure  G.  2(a)  shows  the  case  when  all  four  longitudinal  equations  were  integrated 
using  estimated  values  of  the  lateral  states.  In  Figure  G.  2(b)  the  longitudinal 
values  are  kept  at  the  estimated  values  while  the  lateral  equations  are  integrated. 
In  both  cases  we  see  excellent  agreement  between  the  predicted  and  time  response. 
Figure  G.  2(c)  shows  the  case  when  all  eight  equations  are  integrated  together. 

The  degeneration  in  response  was  dramatic  for  all  the  states.  Another  example 
of  this  sensitivity  is  shown  in  Figure  G.2(d).  This  time  the  response  for  the  same 
maneuver  is  calculated  using  the  wind  tunnel  model  itself.  When  integrating  the 
equations  of  motion,  the  integrated  values  of  a.  and  j8  are  perturbed  by  a  very 
small  deviation  of  0.003  radians  at  254  and  256  seconds.  The  error  introduced 
was  nearly  the  same  as  the  RMS  of  the  measurement  error  on  a  and  /3.  As  can 
be  seen  from  Figure  G.  2(d),  this  very  small  perturbation  produced  extremely 
large  errors  in  all  states.  The  error  is  as  large  as  5°  for  a,  10°  for  j3,  100  deg/ 
sec  for  p,  120°  for  0  and  50°  for  9. 

Output  error  methods  determine  model  parameters  such  that  even  if  the 
models  are  inaccurate  in  some  regions  they  are  compensated  by  models  in  other 
regions  to  match  the  response  for  the  given  maneuver.  From  the  above  examples 
we  see  that  even  the  models  obtained  by  these  methods  can  fail  in  predicting  the 
response  for  a  maneuver  like  the  one  discussed  above  if  it  was  not  used  in  deter¬ 
mining  the  model  since  even  a  very  small  error  in  the  computation  of  a  and  /3 
can  lead  to  a  completely  different  response. 
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6.  PREDICTION  OF  RESPONSES  FOR  A  NEW  MANEUVER 

To  test  how  well  the  model  predicts  an  unknown  maneuver,  control 
inputs  (Figure  6.1)  were  chosen  to  excite  both  the  identified  model  and  the  wind 
tunnel  model.  Sinv soidal  inputs  were  used.  The  amplitudes  of  the  controls 
were  chosen  large  such  that  they  created  a  trajectory  that  did  not  stray 
too  far  into  unmodeled  regions.  From  Figure  6. 1  we  notice  that  the  response 
predicted  by  the  identified  model  and  the  actual  response  are  in  good  agreement. 
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7.  CONCLUSIONS 

The  EBM  methodology  has  been  validated  under  controlled  but  realistic 
simulated  conditions.  The  T-2C  wind  tunnel  model  was  used  to  generate  synthetic 
data.  The  controls,  initial  conditions,  masses  and  thrusts  of  sixteen  actual  T-2C 
flight  test  data  records  were  used  to  excite  the  T-2C  wind  tunnel  model  to  generate 
the  sixteen  synthetic  maneuvers.  The  synthetic  data  were  corrupted  with  realistic 
measurement  noise  levels  and  then  processed  with  an  extended  Kalman-Bucy/ 
Bryson -Frazier  smoother  to  produce  excellent  estimated  values  of  the  states  and 
of  the  forces  and  moments.  Subspace  modeling  together  with  a  step-wise  multiple 
linear  regression  technique  were  used  to  identify  a  global  state/control  dependent 
model  of  the  force  and  moment  coefficients.  The  identified  model  agrees  well  with 
the  T-2C  wind  tunnel  model.  The  high  nonlinearities  were  accurately  identified. 
The  identified  global  model  gave  excellent  predicted  responses  for  a  new  maneuver 
not  used  in  the  identification  process. 
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APPENDIX  A 

EQUATIONS  OF  MOTION  &  MEASUREMENT  EQUATIONS 

A.  1  AIRCRAFT  KINEMATICS 

The  following  kinematical  relations  for  a  rigid  aircraft  are  given  with 

reference  to  body  axes  at  the  center  of  gravity: 

• 

u  =  rv  -  qw  -  g  sin  0+B1 

(A-l) 

* 

v  =  pw  -  ru  +  g  cos  9  sin  0  +  B2 

(A-2) 

w  =  qu  -  pv  +  g  cos  0  cos  ^ 

(A-3) 

p*  -  pq  C41  +  qr  C42  +  q  C43  +  B4  +  B0 

(A-4) 

q  =  pr  C51  +  (r  -  p  )  C52  -  r  C&3  +  B5 

(A  —  5) 

r  =  pq  C61  +  qr  Cg2  +  q  Cgg  +  Bg  +  Cg4  B 4 

(A -6) 

♦ 

0  =  q  cos  <t>  -  r  sin  <j) 

(A- 7) 

♦ 

<b  =  p  +  q  tan  9  sin  ^  +  r  tan  9  cos  $ 

(A -8) 

• 

z  =  -  u  sin  9  +  v  cos  0  sin  i  +  w  cos  0  cos  <j> 

(A-9) 

The  acceleration  terms  B1#  _ ,  Bg  which  are  functions  of  the 

aerodynamic  coefficients  are  given  by  the  expressions: 

^x  pV2 

B,  -  +  -s -  S  C 

1  m  2m  x 

(A-  10a) 

pv2 

B0  =  o -  S  C 

2  2m  y 

(A- 10b) 

pV2  Tz 

B„  =  £ -  SC  +  - 

3  2m  z  m 

(A-  10c) 
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(A-lOd) 


__2  I  I 

f  c,  «rrf£-> 

^  x  x  z  xz 

B5’f-r  Cm  +  fUz  Tx-1.-tm(B3-Tz  )  +  Jz”(Bl-Tx  >!  (A-10e> 
y  ye  — 

B  =£llsb  c  ( _ kjl _ 

6  2  'z  C"  II  -  I  2 

x  z  xz 


m 


) 


The  constant  coefficients  through  Cg4  are  given  by: 
C, 


I  (I  +  I  -  I  ) 
xz  z  x  y 


'41 


'42 


'43 


'44 


'51 


'52 


'53 


I  I  -  I 
X  z  xz 


\  »y  -  y  -  !xz 

■x  ’z  -  4 


I  h 
xz  ex 


I  I  -  I 
X  z  xz 


xz 


xz 


h 


ex 


I  (I  -  I  )  +  r 
_  x  x  y  xz 

'"61  2 - 

I  I  -  I 
X  z  xz 


-  *xz  ^y  "  *z  y 
62  '  ‘ - 


i  i  -  r 

x  z  xz 


I  h^„ 

x  ex 


63  I  I  -  I2 
X  Z  xz 


xz 


'64 


(A-lOf) 


(A-ll) 

(A-12) 

_  (A-13) 

(A- 14) 

(A-15) 

CA-16) 

(A-17) 

(A-18) 

(A-19) 


(A- 20) 


(A-21) 
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A.  2  AIRCRAFT  CONSTANTS  FOR  THE  T-2C 


S  =  23. 7  m 
c  =  2.26  m 
b  =  11.62 


I  =  12212  kg  m2 

I  =  19811  kg  m2 

I  =  25781  kg  m2 

z 


X2  =  0  m 

l  =  .49  m 
z 

e 

The  values  for  the  mass  m, 
given  in  the  following  table. 


h  =  2448  Newtons-m-sec 
ex 

the  thrust  T  and  the  moment  arm  l  are 
x  x 


Table  A-l  Values  of  the  Aircraft  Constants  that  Vary  from 
Maneuver  to  Maneuver 


Maneuver 

Number 

Mass 

(kg) 

Thrust 

(Newtons) 

i 

X 

(Meters) 

1 

5275.3 

4981 

-.006 

2 

5329. 7 

4182 

.031 

3 

5198.2 

3874 

.046 

4 

5039.4 

3683 

.234 

5 

5384.1 

4666 

.018 

;  6 

5039.4 

4240 

.010 

i  7 

4944.2 

4211 

.039 

i  8 

5420.4 

3800 

.047 

9 

5370.5 

3507 

.056 

10 

5325. 2 

3739 

.011 

11 

5216.3 

3658 

.059 

12 

5170. 9 

3458 

.063 

13 

5125.6 

2406 

.212 

1  14 

5080.  2 

2426 

.189 

!  15 

5034. 9 

3707 

.028 

16 

4989.5 

3347 

.052 
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MEASUREMENT  EQUATIONS  FOR  THE  T-2C 

Table  A -2  lists  the  24  variables  that  are  sensed  by  the  T-2C  instrumenta¬ 
tion  package  and  provides  measurement  range  and  location  with  respect  to  the 
normal  flight  test  for  each  sensor. 

The  measurement  equations  relating  the  sensed  variables  to  aircraft 
states  and  controls  are  as  follows: 


-1  /  w  -  qX  +  pyx 
a  =  tan  - ± - —  ) 

m  \  u  +  qza  -  rya  J 


(A- 22) 


_  v  -  p  +  r  x 
8m  tSn  \  u  +  q  -  r 


J-) 

ye  / 


(A -23) 


-  .2 


-  -  v2 


—  .2  .2 


m 


=  [  (u  +  qztr  -  ryj  +  (v  +  rx^-pzj  +  (w  +  pyTr  -  q xj  ] 


(A-24) 


r 


m 


r 


e 

m 


d 


(A- 25) 

(A-26) 

(A -27) 

(A-28) 

(A-29) 

(A-30) 


P 


(A-31) 
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(A-32) 


r 


m 


r 


(A-33) 


<a*  We  '  Bi  '  <r2  4  q2)  *x  +  <pq  '  W  <pr  +  q)zx  (A‘34) 

m  i  i 

where  i  =  NOSE. 

Three  measurement  equations  result  with  k  corresponding  to  c.g.,  NOSE 
and  TAIL  in  the  following  expression; 

<av  *k  =  b2  +  <P9  +  r )  -  <P2  ^  r*2)  Fy  +  (qr  -  P)z  (A-35) 

ym  yk  yk  yk 


Five  measurement  equations  result  with  k  corresponding  to  c.g.,  NOSE, 
TAIL,  R  and  L  in  the  following  expression: 


.  —  *  —  2  2  — 
(az  *k  =  b3  +  <Pr  -  q  )  xz  +  +  p)  yz  -  (p  +  q  )  zz 

m  k  k  k 


5  =  6 

a.  a 

m 


5  =  5 

e  e 

m 


5  =  5 

r  r 

m 


(A- 36) 

(A- 37) 

(A-38) 

(A-39) 


In  Table  A -2,  note  that  x  "4.62,  y  =-.13,  and  z  =0 

XNOSE  XNOSE  XNOSE 

Similarly,  x  =0.12,  y  r  =-.02,  and  z  =0.58. 

c.g.  yc.g.  yc.g. 


.  18. 
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Table  A- 2.  Location  Values  of  Instruments 
in  the  T-2C  Instrumentation  Package; 
Positive  x  means  Instrument  is  Foward  of 
c.g. ,  Positive  y  -  right  of  c.  g. ,  Positive 
z  -  below  c.g. 


Location  in  Meters 
Measured  Quantity  ^  :  = 


Measurement 

Range 


DYNAMICS  RESEARCH  CORPORATION 

SYSTEMS  DIVISION 


Table  A-2.  (Continued) 


Location  in  Meters 

Measurement 

No. 

Measured  Quantity 

X 

y 

Z 

Range 

21 

Normal  Acceleration  at  Left 

Wing  Tip  (a  ) 

Zm  L 

-0.  89 

-5.08 

0.05 

+  3.5  g 

22 

Aileron  (5  ) 

a 

m 

-0.  97 

4.06 

0.05 

+  13  deg 

23 

Elevator  (S  ) 
e 

m 

-5.01 

0.05 

-1.19 

-27. 5  to  15  deg. 

24 

Rudder  (6  ) 

r 

m 

-5.01 

0.05 

-0.72 

+  24  deg. 

24 
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APPENDIX  B 

DESCRIPTION  OF  FILTER /SMOOTHER  ESTIMATOR 

B.  1  DESCRIPTION  OF  THE  STATES 

The  kinematical  relations  (A-l)  through  (A-9)  consist  of  15  states.  These 
are  redefined  as  follows: 


X!=U 

x4  =  p 

x?  =  9 

x2  =  v 

X5=q 

x8=0 

x3  SW 

x6  =  r 

X9  =  Z 

x  =  B 

x.  =  B„ 

x_  „  =  B 

10  1 

16  3 

22 

*13  ’  B2 

*19  =  B4 

*25  '  B, 

The  acceleration  B.  (t),  i  €  (l,  2,  .  .  .  ,  6} ,  is  not  a  quadratic  function 
usually,  but  it  can  be  modeled  by  such  a  function  whose  derivatives  are  allowed 
to  vary  with  time.  In  particular,  it  can  be  modeled  by  a  quadratic  function  whose 
derivatives  are  driven  by  white,  zero-mean  gaussian  noise.  This  type  of  model¬ 
ing  is  adopted  herein.  The  resulting  six  quadratic  stochastic  models  lead  to  the 
following  eighteen  differential  equations: 


X10 

**11  +  ,,10 

I-* 

OO 

“  *14  +  r13 

CO 

=  *17  +  ^16 

X11 

•*12+T,11 

•> 

X14 

*15  +  ^14 

X17 

3  *18  +  *17 

3 

*  0  +  T! 

€ 

=  0  +  T1,  _ 

X,  „ 

II 

o 

+ 

) 

12 

12 

15 

15 

18 

is 
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• 

4 

X19  =  X20  +  r19 

X„„  =  X„„  +  r 

22  23  22 

• 

x„_  =  x„„  +  tu- 

t  25  26  25 

X20  =  X21  +  ^20 

X23  =  X24  +r,23 

X26  =  X27  +  r26 

*21  =  °  +  "21 

*24  =  °  +  "24 

i27  =  0+"27 

where  through  v ^ 

are  white,  zero-mean  gaussian  noise. 

The  above  eighteen  differential  equations  together  with  (A-l)  through 
(A-9)  comprise  the  27-state  dynamic  system  for  which  the  filter  and  smoother 
are  defined. 

Let  x  denote  the  state  vector  (x^  Xg,  ....  x 2?)  and  let  f  (x)  be 
the  non-stochastic  part  of  the  right-hand  sides  of  the  27  differential  equations. 
The  vector  (f  (x)  is  a  column  vector  with  components  f  (x),  f  (x),  .... 

JL  u 

(x).  The  stochastic  dynamic  system  can  be  written  as: 


x  *  f  (x)  +  n  ;  x  (t  )  =  x 

o  o 


(B-l) 


where 


t  »  initial  time 
o 

n  ■  white,  zero-mean  gaussian  process  noise  with 


E  [n  (tj)  n  (tg)]  =  Q  6  (t2  -  tj) 


Q  =  positive  semi-definite  diagonal  matrix  -  contains  the  values  of 
the  magnitude  of  the  white  noises 


x  =  gaussian  initial  condition  vector  with  covariance  P  ;  P  is 
a  positive  semi-definite  matrix  °  ° 


The  vector  n  is  a  column  vector  with  components  r^,  rig,  .  ...  r\^7. 

The  values  of  Q^,  n  -  1,  2,  ...»  9,  are  determined  by  the  computational 
accuracy  of  the  computer.  Section  B.  6  describes  how  the  values  Q^,  n  =  10, 

11,  .  .  .  ,  27,  are  set.  _ _ 


B-2 
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B.  2  DESCRIPTION  OF  THE  OBSERVATION  EQUATIONS 

The  measurement  equations  (A-22)  through  (A-30)  and  (A-34)  through 
(A-36)  consist  of  eighteen  measured  quantities.  The  measured  quantities  in 
(A-31)  through  (A-33)  are  not  considered  herein.  The  eighteen -dimensional 
observation  vector  z  =  (z^,  z^,  .  .  .  ,  zjg)  is  defined  as  follows: 

Acceleration  Group 


*  (a  ) 

z  =  (a  ) 

z_ 

=  (a  ) 

X 

m  NOSE 

4  z 

m  TAIL 

7 

^C.G. 

=  (a  ) 

z  =  (a  ) 

z_ 

=  (a  ) 

z 

m  C.G. 

5  z 

mR 

8 

y 

m  NOSE 

> 

m  NOSE 

Z6=(az  > 
mL 

Z9 

=  ^ay  ) 

TAIL 

Velocity  and  Position  Grout 


z10  '  em 


z12=pm 


z16  *  Vm 


Z11  ^m 


Z13  =  qm 


z17  =  a 

17  m 


Z15  ^m 


z,  .  =  r 
14  m 


Z1fi  = 
lo  m 


Note  that  Equations  (A-34)  through  (A-36)  contain  the  derivatives  p,  q 
and  r.  These  derivatives  are  given  in  terms  of  the  state  x  by  Equations  (A-4), 
(A-5)  and  (A-6).  With  these  substitutions  into  Equations  (A-34)  through  (A-36), 
the  resulting  equations  have  right-hand  sides  that  are  functions  only  of  the  state 
x.  Consequently,  we  can  make  the  following  definitions: 


h^  (x)  *  RHS  of  Equation  (A-34) 

h  (x),  h  (x)  ,  .  .  .  ,  hR  (x)  =  RHS  of  Equation  (A-36)  as  k  varies  from 
J  C.G.  to  L. 

h_  (x),  hQ  (x)  ,  hQ  (x)  =  RHS  of  Equation  (A-35)  as  k  varies  from  C.G.  to 
7  8  9  TAIL 


B 
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where 


The  functions  h1Q  (x),  h.^  (x),  .  .  .  ,  h  g  (x)  are  defined  as  follows: 


(a)  h1Q  (x)  and  h^  (x)  are  the  RHS's  of  Equations  (A-28)  and 
(A-29),  respectively. 

(b)  h  (x),  h^g  (x)  and  h^  (x)  are  the  RHS's  of  Equations  (A-25) 
to  (A-27),  respectively. 

(c)  h^g  (x)  is  the  RHS  of  Equation  (A-30) 

(d)  h17  (x),  h^g  (x)  and  hlg  (x)  are  the  RHS's  of  Equations  (A-22) 
to  (A-24),  respectively. 


The  observation  vector  z  is  measured  at  discrete  times. 


t  <  t  <  .  .  .  <  tXT 
o  1  N 


t^  -  ti_1  =  At  for  all  i  =  1,  2,  ....  N 


At  =  inverse  of  the  data  rate 


where 


The  observation  equation  can  now  be  written  as 


zi~  z  *  k  (x  + 


v  -  gaussian,  zero -mean,  white  noise  sequence  with 


(B-2) 


E  ( v .  v .  )  =  R.  6.  . 
ij  i  1.3 

R^  =  positive  -  definite  matrix  -  covariance  of  the  measurement  noise 


U3 
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SOME  ASSUMPTIONS  AND  DEFINITIONS 

We  assume  that  r  (t),  and  xq  are  uncorrelated,  i.e.,  for  all  t  and  i. 


E  (r  (t)  v  )  =  0 

E  (-  (t)  xT)  =  0 
o 

E  (x  vT)  =  0 
o  l 


We  make  the  following  definitions: 


x.  =  limit  x  (t.  -  c) 


H.  -  (x.) 

1  OX  1 


P.  =  covariance  of  x. 
l  l 


F  (x)  =  (x)  for  all  x 


V.  =  P.  HT 
i  ii 


D.  =  H.  P.  H.  +  R.  (Predicted  measurement  residual  covariance) 

l  ill  i 

=  P^  D^  *  (Kalman  gain) 


For  brevity  of  notation,  we  set 


x.  =  filtered  estimate  of  x  at  time  t. 

l  l 


P.  *  filtered  covariance  of  x  at  time  t. 
i  i 


x.  =  smooth  estimate  of  x  at  time  t. 
i  i 
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B. 4  THE  FILTERING  ALGORITHM 

The  initialization  at  time  t  is 

o 

x  (tJ  = 

o  o 

p  = 
o  o 

The  updates  of  the  state  and  of  the  covariance  at  time  t  is  given  by  the 
Kalman-Bucy  filtering  algorithm: 


4  =  z.  ■  h  (x  J  (predicted  residual) 


x.  =  x.  +  G.  A  z.  (state  vector  update) 
1111 


/\  ip 

P^  =  P.  -  Gi  V.  (optimal  Kalman  update) 


The  propagation  of  the  state  and  the  covariance  from  t.  to  is  given 

by  integrating  the  following  equations: 


x  =  f  (x)  ,  x  (t.)  =  x. 


'«■'  Sf  ,  /^f  ri-'x  "v\T  _ 

P-jjWp  +  tjjS)?)  +  <3 


where 


ssj  ,  .  A. 

p  <y  ■  p4 


The  integrated  values  at  time  t.+^  are  recognized  as  x^  and  p.^; 


that  is. 


Vi =  x  “i+i1 


pi+l =  p  <‘1+1) 


B-6 
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B.  5  THE  SMOOTHING  ALGORITHM 

This  algorithm  uses  the  recursive  form  of  the  discrete  Bryson-Frasier 
smoother,  Bierman  [  6  ]. 

Let  X  denote  the  adjoint  filter  variable  generated  from  a  reverse  time 
Kalman-like  filter.  Define 


X .  =  limit  X  (t.  -  e ) 

1  *  i 

e  -*  o 

X +  =  limit  X  (t.  +  e ) 

l  „  l 

e  -*  0 

The  initialization  of  the  adjoint  variable  X  at  time  t^  is  given  by 


T  -1 

X-t  *  ”  H  D  A  zAT 
N  N  N  —  N 


The  propagation  of  X  from  t.+^  to  L  is  carried  out  by  integrating  the 
following  adjoint  equation  backwards: 


x»-rg(t))x  ,  X  <t.+1>  =  X:+1 


The  integrated  value  of  X  at  time  L  is  recognized  as  X^  ;  that  is. 


X  +  =  X  (t.) 

l  i 


The  smooth  estimate  x.  is  obtained  from  the  equation 


*  A  A  ,  + 

x,  =  x.  -  p.  X . 
i  l  i 


The  update  of  )  at  time  L  is  given  by 

_  A  _  1  rp  P  a 

X  .  *  X  .  -  H.  D.  A  z.  -  H.  G.  X. 
1111  1111 
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B.  6  INPUT  QUANTITIES  FOR  THE  FILTER/SMOOTHER  ALGORITHM 

The  filter /smoother  algorithm  outlined  above  is  a  suboptimal  process. 

It  provides  excellent  results  only  if  the  errors  in  the  estimates  are  small  enough 

so  that  the  predicted  state  at  the  next  observation  time  is  contained  within  the 

linear  region  of  the  true  state.  The  critical  factor  in  maintaining  sufficiently 

small  errors  is  the  setting  of  n  *  10,  11,  ...  ,  27.  The  values  of 

Q  .  nn  =  10,  11,  ....  27  are  set  after  analyzing  the  peak  rates  of  change  of 

B.  (t),  i  *  1,  2,  .  .  .  ,  6.  These  peaks  drive  the  quadratic  model  process.  The 

following  ad  hoc  method  is  used  to  set  the  Q  's:  Let  i  e  {l,  2,  .  .  .  ,  6}. 

nn 

Recall  At  =  1.0  divided  by  the  data  rate.  Define 

« 

A  B.  =  peak  change  in  B^  (t)  over  At  seconds  not  due  to  B^  (t) 

•  «  0  0 

A  B.  =  peak  change  in  B.^  (t)  over  2  At  seconds  not  due  to  B^  (t) 


A  B.  =  peak  change  in  B^  (t)  over  3  At  seconds. 


Let  n  =  10  +  3  (i-1).  We  set 


Q  (At)  =  (A  B  )* 
nn  l 


Qfa+l)<n«)  <2  M  *<‘V 


Q(n+2)  (n+2)  <3  Bi>‘ 


The  computed  values  of  the  Q^'s  that  were  used  in  the  filtering  process 
are  given  in  Table  B-l.  The  values  of  Q^,  k  =  1,  2,  .  .  .  ,  9,  were  set  at 
.000001. 


LbtJlJ 
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Table  B-3.  Covariance  Values  of  the  Initial  States  x^g,  x^,  .  .  .  ,  x ^ 
herein,  n  =  10  +  3  (i-1) 


Quantity 

Values  of  the  Covariance  I 

D 

o 

i  =  1 

i  =  2 

i  *  3 

i  =  4 

i  =  5 

i  =  6 

Po  (n)  (n) 

10. 

.01 

200. 

.01 

4.0 

1.0 

Po  (n+1)  (n+1) 

128. 

4.  2 

500. 

2.0 

55. 

2.0 

Po  (n+2)  (n+2) 

22,  500. 

64. 

10,  000. 

784. 

1200. 

700. 

The  values  of  R  (i.  e. ,  the  covariance  of  the  noise)  are  given  in 
Table  B-4. 


Table  B-4.  Covariance  Values  of  the  Noise  for  the  Measured  Quantities 
Zy  z^.  ...  ,  zlg;  herein,  n  =  1  +  3  (i-1) 


Covariance  Values  of  the  Noise 

Quantity 

§1113 

BB 

i  =  3 

m 

B 

mm 

R 

.0576 

1.1025 

.0576 

.756  x  10"6 

2.8  x  10'6 

.  1849 

nn 

R(n+1)  (n+1) 

1.1025 

1.1025 

.0576 

.756  x  10‘6 

2.8  x  10-6 

7.84  x  10'6 

R(n+2)  (n+2) 

1.1025 

1.1025 

.0576 

2.8  x  10'6 

100 

7.84  x  10~6 

The  best  initial  values  of  the  states  xq  are  obtained  by  averaging  over 
the  trim  portion  of  the  maneuver. 
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STEPWISE  MULTIPLE  LINEAR  REGRESSION  (SMLR) 

SMLR  [  4  ]  is  used  to  obtain  state  and  control  dependent  mode  for 
Cx,  Cy,  etc.  The  SMLR  algorithm  is  presented  briefly  in  this  Appendix 
and  details  can  be  found  in  [4  ].  SMLR  determines  in  stages,  the  best 
value  in  the  least  squares  sense  for  the  coefficients  of  an  equation  of  the  form 
given  below: 

y  =  StVi,Y:*"Vn 

where  y  is  the  dependent  variable 

Xj,  x^ - xn  are  the  dependent  variables 

a  ,  a, , - a  are  the  coefficients  to  be  determined.  For  aircraft 

o  1  n 

identification  y  represents  the  force  or  moment  coefficient. 

THe  independent  variables  are  formed  by  combining  one  or  more  of 

QC  OtC 

the  aircraft  states  and  controls.  An  example  is  a,  0,  ^  ’  2V'  a& 

and  a2. 

The  algorithm  proceeds  with  the  independent  variables  being  added 
one  at  a  time  leading  to  the  following  intermediate  equations: 

y  =  a  +  a.x, 

J  o  11 


y  =  a  +  a  *x  +  a,  'xr 


y  =  ao'  +  ai  "  X1  +  a2  "X2  +  a3 


and  so  on.  Thus  at  each  stage  we  can  obtain  valuable  statistical  information. 
At  each  step,  the  independent  variables  not  used  in  the  regression  are 
inspected  and  the  variable  that  gives  the  largest  reduction  in  the  variation 
of  y  is  chosen.  This  variable  is  then  checked  for  statistical  significance 
by  calculating  the  F-Ratio  of  the  variable  and  comparing  against  the  chosen 
F-Critical  value.  If  it  meets  this  requirement  it  is  included  in  the  solution. 


C-l 
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WIND  TUNNEL  MODEL 

This  section  describes  the  wind  tunnel  model  used  in  generating  the  synthetic 
data.  Wind  tunnel  measurements  of  static  and  dynamic  coefficients  for  the  air¬ 
craft  were  provided  by  NADC  [  5  ].  The  longitudinal  and  lateral  models  will  be 
dicussed  below. 

LONGITUDINAL  MODEL 

The  lift  coefficient  C^,  drag  coefficient  and  the  pitching  moment  coeffi¬ 
cient  C  ^  were  given  as  functions  of  angle  of  attack,  sideslip  angle  and  elevator 

nig 

setting.  These  were  made  available  in  the  form  of  look-up  tables.  For  the 

elevator  settings  of  -25°,  -15°,  -10°,  0°,  10°  and  15?  the  tables  provided 

values  of  CT  ,  CL  and  C  at  angles  of  attack:  -8,  -4,  0,  4,  8,  12,  14,  15,  16,  18, 
L  D  mg 

20,  25,  30,  35,  40,  45,  50,  55,  62,  67,  72,  77  and  83  degrees  while  the  sideslip 

angle  assumed  values  of  0,  2,  5,  10,  15,  20,  25  and  30  degrees  at  every  given 

angle  of  attack.  The  given  model  assumed  symmetry  with  respect  to  angle  of  side- 

o 

slip.  For  angles  of  sideslip  greater  than  30  degrees  the  values  at  8  =  30  were 

used.  The  dynamic  derivatives  C  ,  C  ,  C  and  C  were  given  as 

z  zd  m 

functions  of  only  angle  of  attack.  Tl$3  dependence*  of  these  derivatives  on  the 

angle  of  attack  was  furnished  as  curves.  The  given  model  assumed  negligible 

effect  of  sideslip  angle  on  the  dynamic  derivatives;  and  lateral  controls  on  the 

longitudinal  coefficients.  At  any  angles  of  attack,  sideslip  and  elevator 

deflection,  the  values  of  CT  ,  C_  and  C  were  computed  from  the  tables  using 

L  D  ms 

a  multi-dimensional  interpolation  scheme.  When  the  values  of  a  or  were 
outside  the  given  tables  the  scheme  simply  extrapolated  the  values  from  the 
tables. 

The  axial  force  coefficient  C  ,  the  normal  force  coefficient  C  and 

x  z 

pitching  moment  coefficient  Cm  were  obtained  along  the  body  axes  using  the 
following  equations: 

Subscript  s  denotes  static  value  along  stability  axes. 


P-l] 
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The  lateral  force  and  moment  coefficients  are  first  computed  along  stability 
axes  by  directly  adding  the  contributions  due  to  aileron  and  rudder  deflections. 
These  are  then  transformed  into  body  axes.  The  following  equations  outline  the 
procedure  used  for  this. 


cy  <«.  0.  V  V  =  Cys+  *  ACy(V 


(D.7) 


C*  =  sgn  (1.  3)  C  (  a.  |/9|  .  6  »  0^  =  0) 


(D.8) 


AC  (6  )  =  6  (C  (5  )  -  C  )/12 

y  a'  a  y  a'  ys 


(D.9) 


C  (  a,  3.  5a  =  12.  6^  =  0)  .  5aM;^0 

ys 

■C  (  a.  |j3|  .  =  -12.  6r  -  0) ,  a  0  ;  3  <  0 


C  (6  )  = 
y  a 


Cy  (a,  fi  ,  6a  =  -12.  6^=0),  ^a  <  0  »  0  a  0 

s 


*C  (  a.  |j3|  ,  6&  =  12,  6r  =  0)  ,  6&  <  0  ;  3  <  0 


AC  (5  J  =  16  J  (CJ6)-CJ/25 


y  r 


y  r  y 


CA> 


C  (  a.  3,  6  =  0.  8r  =  -25)  ,  6  <  0;  0  2  0 

ys 

-C  (a,  Ifll  ,  6  =  0.  5  =  25)  ,  6  <  0;  3  <  0 
ya  I  1  a  r  r 


C  (  a.  3  ,  6  =0.  6  =25)  .  6  s  0;  3  2  0 

y  a  r  i 

s 


(  a.  |3  I,  6a  =  0,  6^  =  -25)  >  5  r  a  0;  3  <0 


I-DtJ 


(D.  10) 


(D.  11) 


(D.  12) 


i 
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The  static  coefficients  along  body  axes  are  obtained  from: 
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« 

c  =  c 
y 


(D.  13) 


C„  =  C  „  cos  a  -  C  sin  a. 
In  n 

s  s 


(D.14) 


C  =  C  cos  fl  +  C  a  sin  a 
n  n 

s  s 


(D.  15) 


The  total  lateral  force  and  moment  coefficients  along  body  axes  are  are 
then  given  by: 


*  pb  rb 

Cy  Cy  +  Cyp  zv  +  cyr  2V 


(D.  16) 


c  -  c  *  +  c  P^3  +  C  ^ 
c  1  2V  27 


(D.  17) 


C*  +  C  ^  +  c  SiL 
n  n  2  V  n  2V 

P  r 


(D.  18) 


The  wind  tunnel  model  has  been  compared  against  identified  model  and 
discussed  in  detail  elsewhere.  For  the  sake  of  conciseness  it  will  not  be 
discussed  again  here.  The  salient  feature  of  the  model  was  its  high 
complexity  thus  providing  a  thorough  check  out  of  the  identification  technique 
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Figure  D-3  (Concluded) 


'igure  D-4  Time  Histories  of  Simulated  Maneuvers  (Maneuvers  4) 
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Figure  D-5  Time  Histories  of  Simulated  Maneuvers  (Maneuver  5) 
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Figure  D-10  Time  Histories  of  Simulated  Maneuvers  (Maneuver  10) 
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Figure  D-ll  Time  Histories  of  Simulated  Maneuvers  (Maneuver  11) 
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APPENDIX  E 

ESTIMATION  RESULTS  OF  THE  SYNTHETIC  DATA:  TABLES  AND  CURVES 

Tables  E-l  through  E-3  contain  the  RMS  values  of  the  estimation  errors. 
These  RMS  values  are  obtained  by  taking  the  difference  between  the  noise-free 
measurements  and  the  measurements  obtained  through  substituting  the  estimated 
states  into  the  measurements  equations.  The  RMS  values  in  these  tables  are  for 
Maneuvers  1-3,  5-16.  The  results  for  Maneuver  4  are  contained  in  Section  3. 

The  RMS  residuals  obtained  by  taking  the  difference  between  the  true 
values  of  (t),  B^  (t)  ,  .  .  .  ,  Bg  (t)  and  their  estimates  are  given  in  Tables 
E-4  and  E-5  for  Maneuvers  1-3,  5-16. 

The  comparisons  between  the  true  values  of  B^  (t),  i  =>  1,  2,  .  .  .  ,  6, 
and  their  estimated  values  are  given  in  Figures  E-l  through  E-6  for  Maneuvers 
1  and  3. 
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Figure  E-3  Comparison  Between  the  True  and  the  Estimated  Values  of  Bg  (t) 
and  B„  (t)  for  Maneuver  1 
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Table  E-2  Comparison  Between  Measurement  Noise  and  Estimation  Error  Levels  for 
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Table  F-l  LONGITUDINAL  MODEL  EQUATIONS  FOR  ADDITIONAL 

REGIONS 

1)  -4°  s  a  <  0°,  0°*  |j3|*2°,  -10°  s  6^0° 

C  =  -0.0184 

x 

C  =  0.2371  -  4.856  (a +  4/57.3)  -0.8969  6  -  4.4385  <jc 

2  e  2V 

C  =  0.0302  -  0.9556  (a  +  4/57.3)  -  1. 5100  6 
m  e 

-13.1871  _cjc 

2V 

2)  8°sas  12°,  10"°  s  0  *  15°,  -10°  s  6&  £  -5° 

C  =  0.0562  +  0.9528  (a  -  8/57.3) 
x 

C  =  -0.694  -  5. 5204  (a  -  8  / 57. 3)  -  1.9000  (6  +  5/57.3) 

z  e 

-24.5206 

C  =  -0.0715  -  1.  1750  (a  -  8/57.3)  -  1.8634  (6  +  5/57.3) 

m  e 

-14.6395 

3)  20°  s  a  *  25°  ,  10°  «  |j9  I  *  15°,  -17°  s  6e  ^  "13° 

C  -  -0.047528  -  0.3964  (a  -  20/57.3)  +  0.2553  (  i 3  I  -  10/57.3) 

x 

-0.3815(6  +  13/57.3) 

e 

C  =  -0.08654  -  0.5347  (a -20/57.3)  -  0.7022  (  (01-10/57.3) 

z  _ 

-0.  1086  (6  +  13/57.3)  +  (-2.  2145  +  32.4  a  -  20/57.3) 

C  =  -0.  1070  -  1.478  (a-20/57.  3)  -0.  5959  (  |S | -10/57.3) 
m 

-0.6877  (6  +  13/57.3)  +  (1.2442  -  31.26  a-20/57. 3) 
e 


9S. 

2V 


Table  F-2  Comparison  of  Estimated  and  Wind  Tunnel  Values  of 

.Average  Cv 


a 

0°*  \8\  K  2° 

VI 

o 

CM 

<  5° 

5°£  1/3! 

£  1Q° 

(deg) 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0. 

-0.9218 

-0.7563 

-0.8135 

-0.7238 

t  -0.6387 

-0.8113 

4.0000 

-0.8082 

-0.8336 

-0.8014 

-0.8002 

-0.6387 

-0.8583 

6.0000 

-0.8660 

-0.7993 

-0. 8014 

-0.8002 

-0.6387 

-0.8858 

8. 0000 

-0. 7026 

-0.7649 

-0. 8237 

-0.8537 

■  -0.9093 

-0.9133 

10.0000 

i 

-0.7019 

-0.6518 

-0.9154 

-0.8718 

-0.9030 

-0.9213 

|  12.0000 

1 

-0.6947 

-0.5336 

-0.8991 

-0.8900 

-0.7556 

-0.9293 

|  13.0000 

-0.6643 

-0.6517 

-0.8530 

-0.8566 

-0.7556 

-0,7747 

:  14.0000 

! 

-0.7649 

-0.7648 

-0.7197 

-0.8232 

-0. 4445 

-0.6201 

!  15.0000 

i 

-0. 8437 

-0.6733 

-0.8127 

-0.6381 

-0.5211 

-0.4895 

'  16.0000 

-0.7002 

-0.7104 

-0.7990 

-0.9107 

-0.7517 

-0. 7720 

17.0000 

i 

-0.7316 

-0.7462 

-0. 8012 

-0.8087 

-0.8303 

-0.7441 

18.0000 

-0.9010 

-0.7820 

-0.7759 

-0.7068 

-0.6950 

-0.7162 

19.0000 

-0.9270 

-0. 8565 

-0.7782 

-0.7458 

-0.7409 

-0.6945 

20.0000 

-0.  8468 

-0.9309 

-0.7842 

-0.7849 

-0.6745 

-0.6727 

22.5000 

-0. 9004 

-0.9324 

-0. 6974 

-0.7105 

-0.8154 

-0.7574 

25. 0000 

-0.9783 

-0.9339 

-0. 7037 

-0.6361 

-0.8050 

-0.8422 

27.5000 

-1. 1869 

-0.9998 

-0.7080 

-0.7296 

-0.9233 

-0.8531 

30.0000 

-0.9601 

-1.0656 

-0. 9425 

-0.8231 

-0.0156 

-0.8640 

k 


Table  ^-3  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  C 


a  (deg) 

0°  £ 

0\<2° 

2°  <  \0  !<  5° 

5°<iflj<  10° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0-4 

0. 2653 

0. 1076 

0. 

0.0923 

0.0775 

0.0913 

4-6 

0. 

0. 1046 

0. 

0. 0916 

0.0775 

0.0991 

6  -  8 

0. 1351 

0.  1042 

0. 

0. 0921 

0. 0775 

0. 1033 

8-10 

0.0870 

0. 1038 

0.0782 

0. 0927 

0. 1001 

0.1074 

10  -  12 

0.  1157 

0. 1036 

0. 0814 

0.0930 

0.0481 

0. 1003 

12  -  13 

0. 1148 

0. 1033 

0. 0687 

0.  0933 

0. 

0.0933 

13-14 

0. 1017 

0. 1002 

0. 0662 

0. 0875 

0. 

0.0740 

14  -  15 

0. 0900 

0.  0970 

0.  0727 

0. 0817 

0. 

0.0547 

15  -  16 

0.  0691 

0.0891 

0. 

0.0759 

0.0845 

0.0608 

16  -  17 

0. 0750 

0. 0767 

0.  0594 

0. 0728 

0. 0784 

0.0737 

17  -  18 

0.  0772 

0.0718 

0. 0494 

0.0694 

0.0815 

0.0678 

18  -  19 

0.0795 

0.0668 

0.  0704 

0.0661 

0. 

0.0619 

19  -  20 

0. 1010 

0.0608 

0.0734 

0.  0621 

0.  0407 

0.0611 

20  -  22.5 

0.0705 

0.0548 

0. 0905 

0.  0582 

0. 

0.0602 

22. 5-25 

0. 

0. 0496 

0.  0550 

0. 0537 

0.0482 

0.0545 

25  -  27.5 

0.0354 

0. 0444 

0.0347 

0. 0493 

0. 

0.0488 

27.5  -  30 

0.0205 

0. 0462 

0. 

0.0489 

0. 

0.0459 

30  -  35 

0. 

0.0480 

0. 

0. 0485 

0.0251 

0. 0430 

Table  F-4  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  C 


oc  (deg) 

VI 

o 

o 

i$l<  2° 

Estimated 

Wind 

Tunnel 

0-4 

0. 

0. 1252 

4-6 

0. 

0. 1292 

6-8 

0. 1648 

0. 1318 

8-10 

0. 1209 

0. 1344 

10  -  12 

0. 1388 

0. 1310 

12  -  13 

0. 1237 

0. 1276 

13-14 

0.  1419 

0. 1211 

2°<!/3|<5° 


5°  £  (81  <  10° 


14  -  15 

15  -  16 

16  -  17 

17  -  18 

18  -  19 

19  -  20 

20  -  22. 5 

22.5  -  25 
25  -  27. 5 

27.5  -  30 
30  -  35 


0.  1076 
0.  0829 
0. 

0. 

0.0804 
0. 0916 
0.0534 
0. 

0. 0568 
0. 0361 
0. 


0.  1146 
0.  1077 
0.0935 
0. 0822 
0.0710 
0.  0625 
0.0541 
0.  0544 
0. 0546 
0. 0597 
0.  0649 


0. 

0. 1330 
0. 1330 
0. 1428 
0. 1481 
0. 1350 
0. 1255 
0.  1366 
0.  1518 
0.  1105 
0.  0717 
0.  0613 
0.0312 
0. 

0.  0618 

0. 

0. 0235 

0. 


0. 1440 
0. 1485 
0. 1509 
0. 1533 
0.  1471 
0.  1409 
0.  1342 
0. 1274 
0.  1165 
0.  0946 
0.  0800 
0.0654 
0. 0576 
0. 0498 
0. 0540 
0.0583 
0. 0632 
0.  0682 


0. 1131 
0. 1131 
0. 1131 
0. 1752 
0. 1434 
0.  1498 
0.  1498 
0. 1177 
0. 

0. 

0.  1186 
0.0771 
0. 

0. 0288 

0. 

0. 

0. 0298 
0. 0312 


0. 1563 
0. 1595 
0. 1593 
0. 1591 
0. 1574 
0. 1558 
0.  1601 
0. 1645 
0. 1388 
0.  0966 
0.0827 
0.0687 
0.0592 
0. 0496 
0. 0534 
0.0572 
0.0613 
0.0655 


Table  F-5  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  C 


+ 

a 


a.  (deg) 

0°  < 

\0\<  2° 

2°  <1/5 

!<  5° 

5°  <  8,  <  10° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0-4 

0. 

0.  0462 

0. 

0.  0541 

0. 

0.0538 

4-6 

0. 

0. 0356 

0. 

0.  0505 

0. 

0.0630 

6-8 

0. 

0.0182 

0. 

0. 0263 

0. 

0.  0447 

8-10 

0. 

0. 0008 

0. 

0.0020 

0. 0677 

0.0264 

10  -  12 

0. 

0. 0066 

0. 

0.0272 

0. 

0.0672 

12  -  13 

0. 

0.0123 

0.0812 

0. 0524 

0. 1286 

0.1080 

13-14 

0. 

0.0234 

0. 1382 

0. 0593 

0. 1286 

0.0998 

14  -  15 

0.0713 

0. 0345 

0. 0395 

0. 0662 

0. 0726 

0.0916 

15  -  16 

0.0792 

0.0259 

0. 

-0.0010 

0. 

0.0039 

16  -  17 

0. 

0.0265 

0. 

0. 0023 

0.0674 

0.0272 

17  -  18 

0. 

0. 0300 

0.0974 

0.0220 

0. 

0.0371 

18  -  19 

0. 

0.0335 

0. 

0. 0418 

0. 1663 

0.0470 

19  -  20 

0. 

0.0264 

0. 

0.0394 

0.0737 

0.0347 

20  -  22.5 

0. 

0.0193 

0. 

0. 0370 

0. 

0.0224 

22.5  -  25 

0. 

0.0016 

0.  0700 

0.0174 

0.  0630 

0.0014 

25  -  27. 5 

0. 

-0. 0160 

0. 

-0. 0021 

0. 

-0.0197 

27.5  -  30 

0. 

-0.0284 

0. 

-0.0232 

-0.0524 

-0.0331 

30  -  35 

0. 

-0.0408 

-0. 1041 

-0.0443 

-0.07  28 

-0. 0464 

Table  F-6  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  c 


a  (deg) 

0°<l 

0\<  2° 

2°^i 

t-° 

<  5 

5°  <  jfl[<  10° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0  -  4 

0. 

0. 0467 

0. 2787 

0.0651 

0. 

0.0587 

4-6 

0. 

0. 0315 

0. 

0.0400 

0. 

0. 0335 

6-8 

0. 

0. 0286 

0. 

0.0388 

0. 

0.0332 

8-10 

0. 0410 

0. 0258 

0. 0452 

0.0376 

0. 0397 

0.0329 

10  -  12 

0. 

0. 0226 

0. 0652 

0.0347 

0. 

0.  0277 

12  -  13 

0. 

0. 0193 

0. 

0.  0318 

0. 

0.0225 

13-14 

0. 

0.0179 

0. 

0. 0087 

0. 

0.007  2 

14  -  15 

-0. 0235 

0.  0164 

0. 

-0.  0144 

0. 

-0.0081 

15  -  16 

0. 

0.0370 

0. 

-0. 0346 

0. 

-0. 0186 

16  -  17 

0. 0526 

0. 0257 

0. 

-0.0288 

-0.0570 

-0.0123 

17  -  18 

0.0757 

0.0154 

0. 

-0.0200 

0. 

0. 0005 

18  -  19 

0. 

0. 0052 

0. 

-0. 0112 

0. 

0.0133 

19  -  20 

0. 

0.0033 

0. 0887 

0. 0093 

0. 

0.0305 

20  -  22.  5 

0. 

0. 0014 

0. 1017 

0.  0298 

0. 

0.0476 

22.5  -  25 

0. 

-0. 0022 

0. 

0. 0422 

0. 

0.0581 

25  -  27. 5 

-0. 0514 

-0. 0057 

0. 

0. 0545 

0. 

0.0686 

27.5  -  30 

-0. 1705 

-0.0104 

0. 

0. 0407 

0. 

0.0626 

30  -  35 

-0.  1705 

-0. 0150 

-0. 1506 

0.  0269 

0. 

0.0566 

DYNAMICS  RESEARCH  CORPORATION 

SYSTEMS  DIVISION! 


Table  F-7 


COMPARISON  OF  ESTIMATED  AND  WIND  TUNNEL  VALUES  OF  C 

"r 


Wind 

Tunnel 

Estimated 

a(deg) 

0°  s 

0!<  2° 

2°  ^|/3|<5° 

5°  &|/3  |s  10“ 

0 

0.55 

0. 

5285 

0 

0 

4 

0.53 

0. 

6522 

0.4285 

0 

6 

0.505 

0 

0.4285 

0 

8 

0.48 

0. 

5546 

0. 5334 

0.5908 

10 

0.46 

0. 

3892 

0.4871 

0 

12 

0.44 

0. 

4541 

0.4571 

0 

13 

0.447 

0. 

6364 

0 

0 

14 

0.455 

0. 

3485 

0.4575 

0 

15 

0.450 

0. 

3512 

0 

0.54*1 

16 

0.445 

0. 

5423 

0 

0.4885 

17 

0.39 

0 

0 

0.5738 

18 

0.25 

0. 

2311 

0 

0 

19 

0.275 

0 

0.3687 

0 

20 

0.29 

0. 

2012 

0.2732 

0 

22 

0.3 

0 

0.3287 

0 

25 

0.315 

0 

0 

0 

27 

0.358 

0. 

2912 

0 

0 

32.500 

0.405 

0 

-0.4832 

0 

Table  f-8  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  Ci  + 


a  (deg) 

•°i| 

0\<  2° 

2°sW 

<  5° 

5°  <|sk  10° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0-4 

0. 

0.0155 

0.  0161 

0.  0150 

0.0128 

0.0143 

4-6 

0. 

0.0153 

0. 

0. 0150 

0.0128 

0.0144 

6-8 

0.0174 

0.0149 

0. 

0.0146 

0.0128 

0.0141 

8-10 

0. 0139 

0.0146 

0.0137 

0. 0142 

0.0139 

0.0139 

10  -  12 

0.0142 

0. 0141 

0.0137 

0.0147 

0.  0126 

0.0154 

12  -  13 

0. 

0.0137 

0. 0152 

0. 0153 

0.  0250 

0.0170 

13-14 

0.0083 

0. 0134 

0. 0159 

0.  0162 

0. 0250 

0.0189 

14  -  15 

0.  0127 

0.0132 

0.0179 

0.  0171 

0. 

0.  0209 

15  -  16 

0.0130 

0. 0127 

0.0143 

0.0175 

0.0158 

0.0221 

16  -  17 

0. 0099 

0.0109 

0.0158 

0.0147 

0.  0270 

0.0188 

17  -  18 

0. 0097 

0.  0089 

0. 0154 

0. 0118 

0.0168 

0.0157 

18  -  19 

0.0064 

0.0070 

0. 0093 

0.0089 

0.0143 

0.0127 

19  -  20 

0.0066 

0.0071 

0.0085 

0. 0079 

0.0073 

0.0097 

20  -  22.5 

0.0054 

0.  007  2 

0. 0103 

0.  0069 

0.0079 

0.0067 

22.5  -  25 

0. 

0. 0071 

0.0070 

0. 0075 

0. 

0.0073 

25  -  27.5 

0.0101 

0. 0069 

0. 0058 

0.  0081 

0. 0071 

0.0080 

27.5  -  30 

0. 

0.0038 

0. 0270 

0. 0019 

o. 

-0.0012 

30  -  35 

0. 

0.  0006 

-0. 0167 

-0. 0044 

0. 

-0.0104 

F-9 


Table  F-9  Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  Ci 

6r 


a  (deg) 

0°  <| 

a\<2° 

2°  <  W 

<  5° 

5°  <  jf 

||<  10° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0-4 

0. 

0. 0158 

0. 0148 

0.0167 

0.0143 

0.0178 

4-6 

0. 

0. 0150 

0. 0134 

0.0160 

0. 0143 

0.0180 

6-8 

0.0134 

0.0149 

0.0134 

0.0159  | 

0.0143 

0.0176 

8-10 

0.0138 

0.0147 

0. 0154 

0.0157 

0.0158 

0.0171 

10  -  12 

0.0143 

0.0141 

0.  0147 

0.0143 

0.0161 

0.0143 

12  -  13 

0.0142 

0.0135 

0.  0128 

0.0128 

0. 0068 

0.0114 

13-14 

0.0136 

0.0125 

0.0110 

0.0108 

0.0068 

0.0081 

14  -  15 

0.0176 

0.0116 

0.  0063 

0.0089 

0. 

0.0048 

15  -  16 

0.0152 

0.0100 

0. 0064 

0. 0066 

0. 

0.0016 

16  -  17 

0.0173 

0. 0080 

0. 

0.0061 

0. 

0.0013 

17  -  18 

0. 

0.0074 

0. 0094 

0.0063 

0. 

0.0011 

18  -  19 

0. 0069 

0.0068 

0.0074 

0. 0065 

0.0113 

0.0008 

19  -  20 

0.  0084 

0.  0075 

0. 

0. 0073 

0. 

0.0017 

20  -  22.5 

0.  0082 

0.0082 

0.  0064 

0.  0081 

0. 0035 

0.0027 

22. 5-25 

0. 0052 

0.0063 

0.  01.21 

0. 0065 

0. 0055 

0.0028 

25  -  27.5 

0.0057 

0.0044 

0. 

0.0049 

0. 

0.0030 

27.5  -  30 

0.0034 

0. 0044 

0. 0107 

0. 0080 

-0.0104 

0.0106 

30  -  35 

0. 

0.0043 

II 

0. 

0.0112 

0.0147 

0.0182 

Table  F  10  E  stimated  Values  of  Average  C 'i 


a 

(deg) 

DD 

20^  |S|  <5° 

5°  £  (0|  s  10° 

0-4 

0. 

0. 

0. 

4-  6 

0. 

0. 

0. 

6-8 

0. 

0. 

0. 

8-10 

0. 

0. 

-0.0205 

10  -  12 

0. 

0. 

0. 

.2-43 

0. 

0. 

0. 

13  -  14 

0. 

0. 

0. 

14  -  15 

0. 

0. 

0. 

15  -  16 

0. 

0. 

°. 

16  -  17 

0. 

0.  1315 

0. 

17  -  18 

0. 

0.0576 

0. 1125 

18  -  19 

0. 

0. 

0. 

19  -  20 

0. 

0.0609 

0. 

1 

20  -  22.5 

0. 

0. 

0. 0444 

22  -  25 

0. 

0. 0159 

0.  0270 

25  -  27.5 

0. 

0. 

0. 

27  -  30 

0. 

0.0266 

0. 

I 

30  -  35 


i 


-0.  0446 


-0. 0825 


-0.  1541 
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Table  F-ll 


ESTIMATED  VALUE  OF 


Estimated 

a(deg) 

0°£|j3|<2°  2°  s  j$|<  5° 

5°  *|)3  1  *  10° 

■ 

0 

1.2672 

-2.6910 

0 

0 

0 

0 

0 

1.6017 

1.8868 

6.7775 

4.1287 

6.0588 

5.8826 

10.5690 

10.3835 

24.8933 

23.8968 

14.1890 

11.9654 

0 

0 

-6.6689 

-6.2384 

-9.3653 

-12.6470 

-0.9109 

-1.3264 

-3.9284 

-4.3695 

0 

0 

-2.1175 

0 

-0.4477 

-0.6367 

0.1556 

0.1556 

0.1556 

0 

-  2.3044 
3.6924 
3.6924 
8.6769 
24.7578 
14.5607 
-2.4819 
-7.0379 
-8.9929 
-0.5460 
-5.0747 
0 
0 
0 
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Table  F-12  ESTIMATED  VALUE  OF  C. 


Estimated 


ol  (deg) 

0°  ^jj3|<  2° 

2°  *|0|<  5° 

5°  sj/S^XO  ° 

0 

0 

0.6385 

-0.4331 

4 

0 

0 

-0.4331 

6 

0 

0 

-0.4331 

8 

0 

0 

-0.4331 

10 

0 

0 

0 

12 

14. 1 5t0 

0 

7.2825 

13 

0 

0 

7.2825 

14 

7.7380 

6.2149 

0 

15 

0 

0 

0 

16 

0 

0 

0 

17 

0 

-9.5115 

0 

18 

0 

-8.3432 

0 

19 

-12. 1400 

0 

-19.1830 

20 

-2.4842 

0 

0 

22 

0 

0 

0 

25 

0 

0 

0 

27 

0 

0 

0 

1 _ 30 _ 

0 

0 

2.9820 

Table  F-13Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  C 

n6  + 


cl  (deg) 

VI 

o 

o 

\0\<  2° 

2°<!> 

<5° 

5°<|£ 

t|<  io° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0  - 

4 

-0. 0151 

-0. 0100 

0. 

-0.0097 

0. 

-0.0085 

4  - 

6 

-0.0117 

-0.0088 

-0.0154 

-0. 0086 

0. 

-0.0087 

6  - 

8 

-0. 0047 

-0. 0087 

-0. 0154 

-0. 0080 

0. 

-0.0091 

8  - 

10 

-0. 0085 

-0. 0086 

-0. 0070 

-0.  0070 

-0.  00  91 

-0. 0096 

10  - 

12 

-0. 0092 

-0.0064 

-0. 0078 

-0.0058 

-0. 0101 

-0.0062 

12  - 

13 

-0. 0062 

-0. 0042 

.-0.  0083 

-0.0043 

0. 

-0.0029 

13  - 

14 

0. 

-0. 0021 

0. 

-0. 0026 

0. 

0.0021 

14  - 

15 

0. 0026 

-0. 0001 

0. 

-0. 0008 

0. 0096 

0. 0071 

15  - 

16 

0. 

0. 0010 

0. 

0.  0021 

0. 0078 

0. 0117 

16  - 

17 

0. 

0. 0033 

0. 

0.0063 

0.0128 

0.0152 

17  - 

18 

0. 

0.0035 

0.0089 

0. 0053 

0. 0116 

0. 0125 

1 

00 

rH 

19 

0. 

0. 0037 

0. 

0. 0044 

0.0073 

0.0099 

19  - 

20 

0. 

0.0012 

0. 

0.0002 

0. 0096 

0.0040 

20  - 

22.5 

0. 0048 

-0. 0014 

0. 

-0. 0040 

0.  0040 

-0.0019 

22.5  - 

25 

0. 

-0. 0029 

-0.0080 

-0. 0046 

0. 

-0.  0061 

25  - 

27.5 

-0. 0070 

-0.0044 

-0. 0076 

-0. 0051 

0. 

-0. 0104 

27.5  - 

30 

-0. 1090 

-0. 0013 

-0. 0065 

-0. 0020 

0. 0127 

-0. 0083 

30  - 

35 

0.0067 

0.0018 

0.0065 

0. 0011 

-0. 0076 

-0. 0062 

F-14 


Table  F-14Comparison  of  Estimated  and  Wind  Tunnel  Values  of  Average  Cn 

6a 


a  (deg) 

o°< 

0\<  2° 

9°  <a 

2  <,/3 

<  5° 

s°<  ja|<  io° 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

Estimated 

Wind 

Tunnel 

0-4 

-0.0123 

-0. 0142 

-0. 0192 

-0.0174 

-0. 0166 

-0.0128 

4-6 

-0.0162 

-0.0137 

-0. 0159 

-0. 0168 

-0.0166 

-0.0116 

6  -  8 

-0.0115 

-0.0139 

-0. 0159 

-0.0164 

-0.0166 

-0.0116 

8-10 

-0.0124 

-0.0141 

-0.  0172 

-0. 0160 

-0.0141 

-0.0117 

10  -  12 

-0.0117 

-0.0118 

-0. 0170 

-0. 0119 

-0.0172 

-0.0078 

12  -  13 

-0. 0102 

-0.0096 

-0. 0111 

-0. 0079 

-0.0083  " 

-0.0040 

13-14 

-0. 0076 

-0.0077 

0. 

-0. 0044 

-0. 0083 

-0.0011 

14  -  15 

-0. 0059 

-0.0057 

0. 

-0.0010 

-0.0020 

0.0017 

15  -  16 

0. 

-0. 0055 

0. 

-0. 0030 

0. 

-0.0035 

16  -  17 

0. 

-0. 0038 

0. 

-0. 0038 

-0. 0079 

-0.0076 

17  -  18 

-0.0053 

-0. 0006 

-0.0067 

0.0003 

-0.0124 

-0.0014 

18  -  19 

0. 

0. 0027 

0. 

0. 0044 

0. 

0.  0048 

19  -  20 

0. 

0. 0018 

0. 0076 

0. 0040 

0. 0128 

0.0080 

20  -  22.5 

0. 

0.0009 

0. 0075 

0. 0037 

0. 0163 

0.0112 

22.5  -  25 

0. 

-0.0024 

0. 

-0. 0026 

0. 0090 

0.0055 

25  -  27.5 

-0. 0062 

-0« 0058 

0. 

-0.0090 

0. 

-0.0003 

27.5  -  30 

0. 

-0.0027 

-0.0111 

-0.0083 

0. 

-0.0043 

30  -  35 

-0.0100 

0.0004 

-0. 0222 

-0. 0077 

0. 

-0.0083 

I 
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Table  F-16  ESTIMATED  VALUE  OF  C 

ra 


Estimated 

a  (deg) 

o 

o 

IA 

A 

to 

o 

2°  s:  Ij8  j<  5° 

5°  *ijS|s  10° 

0 

0 

0 

0 

4 

0 

0 

0 

6 

0 

0  i 

0 

8 

0 

0 

0 

10 

0 

0 

0 

12 

0 

0 

0 

13 

0 

o 

0 

14 

2.9110 

0 

-5.4219 

15 

0 

0 

0 

16 

5.2109 

0 

0 

17 

0 

0 

0 

18 

-2.408.' 

-4. 1690 

-6.2156 

20 

0 

0.7845 

1.3860 

22 

0 

0 

0 

25 

0 

-0.8470 

0 

27 

0 

-1.1626 

3.4870 

30 

-1.3187 

-1.8826 

-2.1706 
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Table  F 

-17  LATERAL  PARAMETERS  FOR  ADDITIONAL  REGIONS 

Region  1: 

-4°  *  a  <0°.  0°  *  i0(  *2° 

Region  2: 

8°  *  a  *  12°,  10°  *  |j3|  s  15° 

Region  3: 

20 °  s  a  <  22.5°,  10°  s  \$\*  15° 

Region  4: 

22.5°  s  a  <  25°,  10°  £  Jj3|  s  15° 

Region 


C  C 

?s+  V 

a  a 


0.0007 
-0.1553 
-0.1316  | 
-0.1344  I 


-0.5776 

-0.8150 

-0.8656 

-0.7710  -0.2638 


0.1115 


Region 


C  i  C  !  C  I  C 

y5+  i  y6-j  yp  !  yi 

r  r  i 


0.1617  0.1545 

0.0556 


0.5839  0.6152 

0.6481 


0.1268  -0.0870  0.6106 


Region 
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Table  F-l  7  LATERAL  PARAMETERS  (Cont'd) 


0.0 

-0.0270 

-0.0251 

-0.0246 


-0.1487  -0.1500  -0.1782 
-0.0917  -0.1685  -0.1686 
-0.1860  -0.0794 


-0.1441 


-0. 1084 


Region 


c 

C  . 

V 

V 

r 

r 

0.0230 

0.0161 

- 

-0.0636 

-0.0308 

- 

0.7178  I  0.1567 


-1.2000 

-7.1035 


3.3587 
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Table  F-17  LATERAL  PARAMETERS  (Cont'd) 


Region 

C 

n 

0 

C 

na 

c 

* 3 

1 

o 

• 

o 

- 

0. 1098 

2 

0.0160 

- 

0.0130 

3 

-0.0007 

0.0281 

-0.0105 

4 

0.0 

0.0351 

0.0 

Region 

Cn6  + 
r 

V 

r 

c 

n 

P 

1 

-0.0383 

-0.0593 

- 

2 

-0.0600 

-0.0575 

-a.  0663 

3 

-0.0232 

i -0.0199 

-0.0577 

4 

-0.0298 

; -0. 0228 

-0.0652 

0.0182  -0.0090 

0.0105 


2.6847  i 
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APPENDIX  G 

This  appendix  contains  the  predicted  responses  and  the  comparisons 
with  tiie  wind  tunnel  responses  for  Maneuvers  1-3,  5-10,  12-16  and  the  first 
half  of  Maneuver  11.  Note  that  four  pages  of  plots  are  given  for  each  maneuver. 
Each  page  contains  two  states.  Altogether  the  responses  for  the  states 
a,  /3,  V,  p,  q,  r,  9  and  <t>  are  given.  The  solid  line  is  the  response  of  the 
wind  tunnel  model  and  dotted  lines  with  crosses  represent  the  response  of  the 
identified  model.  All  the  response  plots  are  given  in  Figure  G.l.  Figs.  G.  2(a) 
throughG.  2(d)  demonstrate  the  sensitivity  of  the  aircraft  model  response  to 
small  deviations  in  the  states . 


Ti 


TRUE  VALUE 


PREDICTED  VALUE 


.00  12.00  IS. 00  20.00  24.00 

TIME 


TRUE  VALUE 


PREDICTED  VALUE 


TIME 


T I ME  C SEC ) 


Figure  G.l  Continued  (Maneuver  8) 
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Figure  G.l  Continued  (Maneuver  10) 
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